Spin-orbit interaction in the k • p theory for cubic crystals 
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Recent work by Elder, Ward and Zhang [Phys. Rev. B 83. 165210 (201 1)] has shown need for correction and 
modification of current implementation of the k • p method and operator ordering scheme using the interaction 
parameters defined under double group consideration. This manuscript examines the difference in treatment of 
spin-orbit interaction under the single and double group formulations. We show that the restriction to the adapted 
double group bases, brought about by the imposition of single group selection rule in calculating the k ■ -k 
interaction, is not appropriate. In addition, the unitary transformation employed in the literature to diagonalise 
the intra-band spin orbit interaction in semiconductors with diamond lattice can not remove the inter-band terms. 
It leads to a bases set for valence band ordered differently from D^ in the 0(3) group thus invalidating any 

2 

correlation of magnetic quantum number to the z component of angular momentum. Under the double group 

consideration, spin-orbit interaction affects all the zone centre states and offers a mechanism for changing the 

relative positions of various zone centre states in the conduction band. In addition to the mixing caused by k 

independent spin orbit terms, formation of hybridised orbitals under double group rules also produce mixing. 

This can lead to inversion in materials such as a-tin. The re- arrangement of conduction band zone centre states 

leads to a negative 72 in most materials with respect to the valence band bases having the same order as D^ in the 

2 
0(3) group. The double group formulation is also required in the description of Zeeman interaction, including 

the S ■ B term, for the mixed zone centre states. It provides a correspondence between second order interaction 

parameters and the Luttinger invariants. A form of the Hamiltonian generally used in the literature is obtained 

from the results of EWZ using a unitary transformation of the valence band basis and re-definition of some of 

the Luttinger invariants. Such a unitary transformation enables explanation of optical transition selection rules 

for this form of Hamiltonian and has the effect of inverting the signs of 72 in the valence band. The inter-band 

block of the transform Hamiltonian differs from that of EWZ but the associated change of sign of 72 ensures the 

same coupling. Various single group formulations are compared to the double group formulation with the the 

critical difference identified as mixing caused by inter-band spin orbit interaction. The double group formulation 

is used to model the dispersion in Ge and Si. It overcomes deficiencies in single group formulation and is able 

to explain the conduction band/spin split off band effective mass concurrently using one set of parameters. 

PACS numbers: 71.7().Ej, 71.15.-m, 71.20.Mq, 71.20.Nr, 71.70.Fk 



I. INTRODUCTION 

The k • p method has been developed under the single 
group formulation, initially without consideration of the spin 
degree of freedom' ^ (referred to as the DKK model after 
Dresselhaus-Kip-Kittel^ in this manuscript). The effects of 
electron spin has been considered through the spin orbit inter- 
action terms in the one electron Hamiltonian, obtained from 
Dirac equation using Foldy-Wouthuysen transform^. The 
treatment of spin orbit interaction differentiates the imple- 
mentation of fe • p method into single and double group for- 
mulations. In this manuscript, a single group formulation of 
the k ■ p method utilises single group selection rules and the 
group theoretical method in the evaluation of k ■ p/k ■ n per- 
turbation, whereas the double group formulation make use of 
double group selection rules in the evaluation of the k ■ tt per- 
turbation. In constructing models with a limited near set of 
zone centre bases, the method of Lowdin"* is frequently used 
to account for the effects of remote set of states as Lowdin in- 
teractions. This method can be derived from quasi-degenerate 
perturbation theorjEEl 

To incorporate the effects of electron spin, double group 
selection rules are required, as energy eigenfunctions of the 
unperturbed Hamiltonian form the bases of irreducible rep- 
resentations (IR) of the double group. Recent work by Elder, 



Ward, and Zhang^ (referred to as EWZ in this manuscript) fol- 
lows this prescription, both in terms of classification of zone 
centre energy eigenstates, and the use of selection rules in the 
evaluation of the matrix representation of the k-n perturbation 
with respect to the zone centre states. The use of double group 
selection rules enables the work of EWZ to take into account 
mixing produced by the inter-band k independent spin orbit 
interaction under the one electron theory. From the perspec- 
tive of a semi-empirical technique, it also accounts for mixing 
due to many electron effects such as hybridisation in the for- 
mation of covalent bond in cubic semiconductor crystals. 

The method of invariant featured prominently in the de- 
velopment of the k ■ p method** '^. Without specifying the 
bases, the Hamiltonian is obtained by ensuring its invariance 
under symmetry operations of the lattice. It is expressed in 
terms irreducible perturbation components such as wavevec- 
tor k, corresponding generators, and a set of invariants com- 
monly known as Luttinger invariants. This analysis can be 
performed using either single or double group selection rules. 
Correspondence has been established between the DKK, and 
EWZ models derived from perturbation theory and models 
derived under the method invariant using single, and double 
group selection rules, respectively. Under double group se- 
lection rules, there are five invariants describing the valence 
band, two for the spin orbit split off band and three for the 
inter-band block between the valence band and spin split off 



band. In contrast, there are only four invariants describing the 
corresponding zone centre states under the single group selec- 
tion rules. These invariants are expected to be linearly inde- 
pendent of one another and relations formed between invari- 
ants and second order interaction parameters obtained under 
perturbation theory should preserve such linear independence. 

Historically, spin-orbit interaction is considered as a pertur- 
bation under the single group formulatiorPES Various mod- 
els have been developed to account for experimentally ob- 
served effects of spin. These include the relativistically cor- 
rected DKK model** '", the Kane model'^ '■*, extended Kane 
model'^ '*, and the Weiler modeP^. They differ in the classi- 
fication of the zone centre states used in the perturbative ex- 
pansion, but all make use of the single group selection rules in 
the evaluation of fc • p or fc • tt perturbations. Implicit in these 
models is the absence of mixing between zone centre states 
without spin during the process of constructing zone centre 
states with spin. The complexity of these models increases 
towards that of the double group formulation, with the Weiler 
model being able to explain the finite invariant q and gener- 
ally different Luttinger parameters in blocks coupled by the 
spin orbit interaction. However, all these models overestimate 
the zone centre conduction band effective mass in Ge, and 
underestimate the effective mass of spin-split off band in Si. 
There are similar issues in compound semiconductors where 
Hermann-Weisbuch"* parameters are introduced to attribute 
such discrepancies to unknown remote states. It is generally 
assumed that the method of invariant and perturbation theory 
approach refer to the same zone centre bases set when com- 
paring them. 

There are three elements in the treatment of spin-orbit in- 
teraction under the framework of single group formulation; 1) 
there is an assumption that the spin-orbit interaction is small 
allowing, a) the use of perturbation theory, b) the use of single 
group product bases as 'zero' order un-perturbed wave func- 
tions, and c) the neglect of the k independent inter-band spin 
orbit interaction which causes mixing; 2) the spin orbit in- 
teraction matrix for the valence band referring to the product 
bases is obtained from 'p' like atomic orbitals; 3) the unitary 
transformation for the degenerate valence band, which diago- 
nalises the intra-band k independent spin-orbit interaction, is 
also derived from the 'p' like atomic orbital. This manuscript 
focuses on each of these elements and shows that some of 
these assumptions are inappropriate for the consideration of 
spin-orbit interaction of electrons in cubic crystals. It is found 
that the single group selection rules imposes restrictions upon 
the bases. The existence of a finite spin orbit splitting is well 
known experimentally, and requires a finite intra-band k inde- 
pendent spin orbit interaction when evaluated against adapted 
double group bases of the single group formulation. This im- 
plies that the inter-band k independent spin orbit terms should 
also be finite, and poses a contradiction between the mixing 
caused by such terms, and the restriction placed upon adapted 
double group bases from adherence to the single group selec- 
tion rules. 

This manuscript is divided into twelve sections. Section|ll[ 
III examines the treatment of spin orbit interaction under cur- 
rent single group formulation, and discusses the limitation of 



the adapted double group bases when used in the single group 
formulation. Section |IV] examines the consequence of a uni- 
tary transformation which permits the derivation of the 6-band 
Hamiltonian widely used in the literature. The group theoreti- 
cal background of the unitary transformation, both in terms of 
choice of representation matrices and order of basis, are dis- 
cussed in appendixlFl SectionlVldiscusses the different mech- 
anisms of mixing arising from the construction of zone centre 
cell periodic wave functions, that subsequently require the use 
of double group formulation. This further illuminates how the 
spin orbit interaction affects the ordering of conduction band 
zone centre states derived from anti-bonding hybridised or- 
bitals. Sections Vll|VII introduce the treatment of magnetic 
interactions under the double group formulation. The corre- 
spondence between the Luttinger invariants from both sym- 
metric and anti-symmetric Lowdin terms and second order in- 
teraction parameters in each block are established. The de- 
pendence of the Hamiltonian on the orientation of magnetic 
field is also derived. The following three sections then make 
a detailed examination of the treatment of spin orbit interac- 
tion under the single group formulation with the use of the 
group theoretical method. Section VIII describes the treat- 



ment of spin orbit effect under the DKK model using relativis- 
tic corrections. The matrix representation of H^^i and -ff|o2 
are related to the matrix representation of the p operator. Sec- 
tion IX extends the Weiler model to incorporate remote states 
of all symmetries. Comparisons are then made between the 
Weiler model and double group formulation of EWZ, with 
the main distinction being the allowed mixing of single group 
states under EWZ. Section |X] re-examines the Lowdin inter- 
action obtained from quasi-degenerate perturbation theory. It 
shows that the double group bases, used in EWZ, are the log- 
ical choice in presence of inter-band fe independent spin orbit 
interaction. Any form of Hamiltonian obtained from method 
of invariant refers to zone centre bases of the double group 



formulation. Section XI extract double group interaction pa- 
rameters for Si and Ge from experimentally derived Luttinger 
parameters, and compares dispersion relations obtained near 
the zone centre under the different models. The last section 
summarises the differences between the double group formu- 
lation and other models in the literature. The key difference 
between the approaches is mixing of zone centre states caused 
by k independent spin orbit interaction, and consideration of 
hybridised orbitals with spin degree of freedom. Other gen- 
eral points that have arisen from the development of double 
group formulation are also discussed. 

In this manuscript, the notation of Koster et al.CSI on group 
representations are used. The Condon Shortley phase con- 
vention is followed when spherical harmonics or irreducible 
spherical tensor operators are used. Most of the discussions 
focus on materials with a diamond lattice, though some re- 
marks are made in relation to materials with a zincblende lat- 
tice. To facilitate discussions, the first order interaction be- 
tween states belonging to various IRs are given in appendix 
[B] together with the bases used in appendix IA] The DKK 
Hamiltonian and relations between Luttinger invariants due 
to Lowdin term and second order interaction parameters, are 
obtained under the single group formulation in appendix ID] 



Similar relations between the Luttinger invariants and second 
order interaction parameters under the double group formula- 
tion are also obtained using perturbation theory for the sym- 
metric Lowdin interaction, and shown in appendix ID] The re- 
lationship between some of the first order interaction parame- 
ters under single group formulation are obtained in appendix 
lElby transforming the first order interaction with respect to 
the single group product bases. A distinction is made be- 
tween Luttinger invariants obtained from method of invariant 
and experimentally determined Luttinger parameter. The lat- 
ter is dependent on the model used to obtain the parameters 
from experimental observations. 

A remark is also appropriate here on the terms used in the 
description of zone centre basis functions for the perturbative 
expansion, and as bases of IR of the symmetry group of the 
unperturbed Hamiltonian. The solution of the DKK Hamilto- 
nian at fc = forms the basis of single group IRs. Conversely, 
the transformation properties of eigenfunctions of the unper- 
turbed DKK Hamiltonian, may be represented by the relevant 
single group basis functions of the IR. The single group prod- 
uct bases, obtained by direct product of spinor states and bases 
of one single group IR, form the bases of double group rep- 
resentation which is generally reducible. The unitary trans- 
formation of a single group product basis into adapted dou- 
ble group basis decompose such representation into IRs of the 
double group. It is able to represent the transformation prop- 
erties of energy eigenfunctions of the unperturbed Hamilto- 
nian, which does not contain mixed states derived from single 
group states with different energy. In other words, the unper- 
turbed Hamiltonian does not contain k independent inter-band 



interactions. The use of single group product bases or adapted 
double group bases does not imply treatment of spin orbit in- 
teraction unless the relativistic corrections are included. The 
double group bases can represent the symmetry properties 
of eigenfunctions of Hamiltonian containing k independent 
inter-band interactions. Due to multiplicity in the decompo- 
sition of direct product of double group IR Fg , two set of 
independent bases in the Fg IR are required to fully describe 
the angular dependence of relevant interactions/operators. For 
this purpose, the adapted double group bases from the same 
single group IR are not independent. The spin operator S 
works on the spin degree of freedom and has matrix represen- 
tation with respect to the zone centre bases. With respect to 
the single group product bases, the matrix representation of 



S is simply |ct 
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where a^ are Pauli matrices, and 



n is the dimension of the single group IR. The matrix repre- 
sentation of S with respect to the adapted double group bases 
may then be obtained using the relevant unitary transforma- 
tion. This is, however, not the case when the matrix represen- 
tation of S is formed between general double group bases. It 
can be obtained using group theoretical technique by consid- 
ering the transformation properties of S. 
II. SPIN-ORBIT INTERACTION UNDER THE SINGLE 
GROUP FORMULATION 



The time independent one electron Schrodinger equation in 
absence of magnetic field may be cast into the following form 
after application of the Bloch's theorem. 
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The zone centre states for the perturbative expansion are ob- 
tained from the solution of the Hamiltonian i/o, by se tting 
k — ko. This problem can then be treated in two waygSCSl 
The representation is formed at the point of expansion (fco = 
in this case) with H^^^ incorporated as part of the unper- 
turbed Hamiltonian Hq, which yields the double group formu- 
lation of EWZ. This requires the use of double group selection 
rules. The zone centre bases are therefore the eigenstates of 
Hq, and form the bases of IR of the double group. In the 
second approach, the product between the solution of Hq and 
spinor state are used as bases, but the representation of H^^i 
is not diagonal, and both intra-band and inter-band terms ex- 
ist. An additional assumption is then made that single group 
selection rules are applicable when treating the k • n inter- 
action. From group theoretical viewpoint, the two approaches 



differ in the use of single or double group selection rules when 
dealing with the k • n perturbation. They are formally equiv- 
alent under one electron theory, and should yield the same 
dispersion relation, if the bases for perturbation expansion are 
complete. The second approach thus assumes single group se- 
lection rules (for the k ■ tv perturbation), requires that the one 
electron Schrodinger equation is valid, and give rise to the 
single group formulation generally used in the literature. The 
advantage of using the second approach is the link between 
Luttinger invariants among bands which are degenerate under 
DKK'-*. 

In constructing the Hamiltonian with a limited set of bases, 
a unitary transformation of the single group product bases 
is made to remove the weak interaction between the near 
set and remote set to the desired order in k (see section 



[X| . The unitary transformation should diagonalise the unper- 
turbed Hamiltonian Hq which contains H^^i. However, the 
single group selection rule imposes the requirement of no 
'mixing' between different non-degenerate zone centre states 
of the DKK model^'^. Therefore, the unitary transformation 
occurs within the space spanned by each single group prod- 
uct bases, and it preserves the orthogonality between states 
derived from non-degenerate solutions of Hq. In absence of 
the inter-band k independent spin orbit interaction, the unitary 
transformation which diagonalises H^^i also diagonalises Hq. 
Under these circumstances, the adapted double group bases 
under the restrictions of the single group formulation, can rep- 
resent the angular dependence of the zone centre states. How- 
ever, in the presence of inter-band k independent spin orbit 
interaction, there is a conflict between the use of single group 
selection rules and the need to diagonalise H^^i . This follows 
from the fact that unitary transformation, which occurs within 
the space spanned by single group product bases, can not di- 
agonalise the inter-band k independent spin orbit interaction. 
The validity of single group formulation within the framework 
of one electron theory then rests on whether the inter-band k 
independent spin orbit interaction can be neglected. 

Under the DKK model with relativistic corrections, the 
Hg^i term is treated as a perturbation in parallel with the k ■ p 
perturbation. Using the single group product bases to repre- 
sent zone centre states, this perturbation can be separated into 
the spatial and spin parts such that. 
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where, Sa.p — f c and £,„„ are the matrix representation 
of S and WV{r) x p eval uated with respect to the spinor 
states and single group basi ir^ l ^ , respectively. The operator 
\/V{r) X p transforms according to FJ IR under the Oh group 
symmetry operations. The matrix representation of the opera- 
tor -ff|oi with respect to the single group product bases can be 
constructed from a tensor product of Sq/3 and£,„„. Selection 
rules derived from group theory are applied when evaluating 
this matrix using Wigner-Eckart theorem. They indicate that 
spin-orbit interaction can only affect zone centre states asso- 
ciated with single group states with r4 and Fg symmetry as 
'intra-band' terms, and states with F^,FJ, and F^ symme- 
try as 'inter-band' terms. For zone centre states with T'^ and 
F4 symmetry in the valence and conduction bands, the uni- 
tary transformation with single group restriction diagonalises 
the intra-band terms, and partitions the representation space as 



-r^^ 



. TJ and F2 



'Tr^Fg- 



) Fg . It is impor- 



tant to note that while the intra-band k independent terms are 
diagonalised, the inter-band terms remain off-diagonal. Such 
inter-band terms, particularly those which mix single group 
states with different symmetry, contribute to Lowdin terms 
via their modification of zone centre states. These particular 
terms are, however, neglected in all single group formulations 
of fc • p theory. 

The angular dependence of the representation matrices of 
H^^i evaluated with respect to these bases are shown in Ta- 
ble|l] The unitary transformation which takes the single group 
product bases into the adapted double group bases for the 
F^ ® r^(Fj|") or Fg ® r(r(F4 ) IRs, are listed in appendix 



[C] They are obtained from relations between the product basis 
and the corresponding double group basis published by On- 
odera et al.l^S or those given in EWZ. Both similarity trans- 
formations can diagonalise the k independent spin orbit in- 
teractions evaluated with respect to the single group product 
bases that involve F4 and F^ IRs, as shown in Table [ll 

There are two issues related to the matrix representation of 
H§gi and the unitary transformation under the single group 
formulation. First is the reduced tensor elements of the matrix 
representation of H^^i. For the states with F^ and F4 sym- 
metry in the valence and conduction band, the reduced tensor 
element of the intra-band term is generally equated to A/3, 
where A is the relevant experimentally observed spin splitting 
between the related states under double group classificatiorP^l 
The spin splitting is clearly finite and the intra-band term then 



can not be neglected. This is examined further in section VIII 



where the matrix representation of H^^i is related to matrix 
representation of momentum in the same way for both inter 
and intra-band terms. If one accepts a finite reduced tensor el- 
ement for the intra-band term, one also expect that the reduced 
tensor elements of the symmetry permitted inter-band terms 
should also be finite. Therefore, inter-band k independent 
spin orbit interaction may not be neglected. This raises the 
important issue of what the differences are between the prop- 
erties of zone centre bases of the double group formulation, 
and those of the adapted double group bases under the restric- 
tion of the single group formulation. As we shall see in the 
next section, the adapted double group bases with the restric- 
tion of the single group formulation are incomplete. Mixing 
due to inter-band k independent spin orbit interaction is cru- 
cial to the difference between the double group formulation 
and single group formulation under the one electron theory. 
Many electron effects to be discussed in section |V] also lead 
to mixing making the use of the double group formulation a 
necessity. 

The second issue relates to the order of bases within the 
double group IR. Historically, the valence band and spin orbit 

split off band zone centre states are considered to form the ps 

2 

and pi mutliplets^ even though the lowest basis of Fj^ IR has 
I — 2. The same unitary transformation is applied to states 
with F4^ and T^ symmetry. Consequently, the historically 
chosen transformation takes the single group product basis 
from F4 into basis of Fg" Fg (F4^) according to the | J, AIj) 
notation. When the same transformation is applied to the sin- 
gle group product bases involving V^ IR, such as zone cen- 
tre states of valence band of Ge, the resulting adapted double 
group basis clearly form representation of F^ F^(Fj^) but 
is not ordered according to the \J, Mj) notation. This means 
that it is inappropriate to interpre the AIj quantum number in 
terms of z-component of angular momentum. The impact of 



this discrepancy is discussed further in section IV 



A final observation is the energetic ordering in the valence 
band due to spin orbit interaction, obtained from the specific 
lowest single group basis and shown in Table IT] is in contra- 
diction to experimental observation. In fact, the sign of the 
reduced tensor element is dependent on the / number of the 
single group bases based on the central field approximation. 
So the results shown in Table 11] remains compliant to group 



TABLE I. Spin orbit perturbation calculated using Eq.([2]l. 
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theoretical method. 



III. INCOMPLETENESS OF ADAPTED DOUBLE GROUP 
BASES UNDER THE SINGLE GROUP FORMULATION 

In this section, the use of adapted double group bases within 
the single group formulation is examined to establish whether 
such basis are complete. The adapted double group basis, 
obtained by unitary transformation of single group product 
basis, clearly form the basis of representations of the dou- 
ble group. The matrix representation of a given operator, 
with well defined transformation properties, can be evaluated 
with respect to the bases of the IR. The number of linearly 
independent matrices required to describe such an operator 
is determined by the selection ruleP^HlIl Multiplicity in the 
decomposition of the direct product of Fg IRs of the Oh 
group is a key property that differentiates double group IR 
from the single group IR!22J|26| Pqj- double group IR of F^, 
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while in contrast for single group IR of Fg it is, ± g w ^ 5 
{rf © F;[ ® F^} ® [F^jlSl. According to appropriate se- 
lection rules, operators such as total angular momentum J, 
which transforms according to the FJ IR, are described by 
two linearly independent matrices when evaluated with re- 
spect to the double group bases of the F^ IR, but only one 
linearly independent matrix when evaluated with respect to 
the single group bases of F^ IR. For the total angular mo- 
mentum operator J, one must be able to obtain two linearly 
independent matrices by means of enumeration of different 
bases of the T^ IR. Thus, the double group bases are only 
complete if multiple set of bases yield the required number of 
linearly independent matrices. 

Since the Oh group is a sub group of 0(3), the bases of 
IR of the 0(3) group also form bases of representation of 
the Oh group and are generally reducible. Similarly, sin- 
gle group bases may be constructed from bases of IR of the 



50(3) grouf)2SI22l The adapted double group bases can then 
be constructed from the bases of single group IR using the 
transformations described in appendix O Letting V";'^ define 
the relation between the bases 1 0^) of single group IR F, and 
bases li/;') of the Di IR of 50(3) group, we have. 
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where I is the orbital angular momentum quantum number 
Then the adapted double group bases may be related to the 
product of spin states with bases I?/;') of the Di IR of the 
50(3) group, by a matrix defined by. 



U[ = (I2 <» vf) U^ 



(3) 



where U^ is the unitary transformation which block diago- 
nilses the representation matrices, and are given in appendix 
[C] The matrix representation of the angular momentum oper- 
ator evaluated with respect to the adapted double group bases 
can be obtained from. 
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where J'^'p^ is the matrix representation of J operator eval- 
uated with respect to the product bases between the spinor 
states and bases of Di IR of the 50(3) group. The matrix rep- 
resentation J'sp'' (in unit of K) can be readily obtained from 



Pauli matrices and L 
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which is the matrix rep- 



resentation of the orbital angular momentum evaluated with 
respect to bases of D 
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l(2;+l)x(2; + l) + l2x2 ® i(2i + l)x(2/+l)- (5) 



Considering the lowest basis of Fg IR with / = 2, the 
Cartesian components of the matrix representation for the to- 
tal angular momentum (in units of h) in the F^ partition can 
be obtained using Eq.lH]), and are given as. 
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These matrices constitute one of the two linearly independent 
sets of matrices required to describe the angular dependent 
part of the matrix representation of the total angular momen- 
tum operator. If the adapted double group bases from one sin- 
gle group IR (r^ in this case) are complete, the second set of 
linearly independent matrices should be obtainable from other 
bases with different I number. However, results from using 
single group bases with I = 4 (see Table |IX| l or any other 
bases transforming according to F^ IR, returns a set of ma- 
trices linearly dependent on those given in Eq.dSll. This is not 
surprising given the way in which the adapted double group 
bases are constructed from bases of one single group IR. The 
single group selection rules stipulate that there is only one lin- 
early independent matrix in the representation of J with re- 
spect to the single group bases. The direct product of spinor 
states with such bases followed by the unitary transformation 
can not generate any linear independence among the matrix 
representation of J with respect to single group product bases. 
One may also note that the matrix representation J'' shown 
in Eq.lpb, does not satisfy the condition exp[— ij7'47r] = lEI] 
This reflects the fact the adapted double group bases of F^ and 
F^ remain coupled. Since zone centre states form bases of IR 



of the double group, J'' must be scaled accordingly to satisfy 
the condition exp[— Zj/Mtt] = 1. Under the single group for- 
mulation, the matrix representation of the angular momentum 
operator evaluated with respect to the adapted double group 
bases of the F^ valence band may always be represented by. 
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■jadapt^^-^ = 1 or Ci e Z. 



where J'^ = — SJ"^ and exp 

Therefore the adapted double group bases with restriction to 
one single group IR parentage are incomplete due to their in- 
ability to generate the required number of linearly indepen- 
dent matrices through enumeration. 

The representation of FJ may also be generated from the 
direct product F+ (g) F+ = F+ or F+ (g) F+ = F+ ® F+. 
The second set of linearly independent matrices can be eas- 
ily generated using either the general double group bases with 
well defined j number, or adapted double group bases from 
different single group parentage. To illustrate this point, the 
angular dependent part of matrix representation of J^ is eval- 
uated with respect to the basis of Dt IR of the 0(3), which 

2 

also form basis of F^ IR of the Oh group, are given by. 
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Alternatively, they may be obtained from adapted double 
group bases derived from single group F^ IR. This set of ma- 
trices are linearly independent from those obtained in Eq.(J6|, 
but both are obtained from basis that transform according the 
F^ IR. The matrix representation of the total angular mo- 
mentum operator, evaluated with respect to the general double 
group bases under the double group formulation, is given by. 



J, 



double 



ClJu + C2 Ju- 



(9) 



The 

exp 



coefficient ci and C2 are subject to the condition 
— z J!f''"^''^47r) = 1. This requirement yield 

1 

8 
1 

8 
It is clear that the ratio of the eigenvalue 

|, I) and ||, i) states is (3c2 — Ci)/(c2 — 



ci 



C2 



(3n — m) 
(3m — n) 



(10a) 



(10b) 



where m,n € J. 

of Jdouble foj. jjjg 



3ci) = m/n for the general double group bases but a constant 
1/3 for the adapted double group basis, as indicated in Eq.(J7|. 
This particular feature of the zone centre energy eigenfunc- 
tions under the double group formulation, namely requiring 
two linearly independent matrices to represent the angular de- 
pendence of J, can never be represented by the product bases 
from one single group IR or the corresponding adapted dou- 
ble group bases. Thus the difference between the single, and 
double group formulations can be attributed to mixing of dif- 
ferent non-degenerate single group states in the construction 
of compatible zone centre states under the double group clas- 
sification. 

The matrix representation of j^ with respect the zone cen- 
tre states, as bases of F^ IR, is also diagonal. Hence, the zone 
centre energy eigenstates are always eigenstates of total angu- 
lar momentum, but the eigenvalues of J^ and J^ are not gen- 
erally integer multiples of /i^/4 and h/2 respectively. They 



are material dependent, and the eigenvalue of j^ is propor- always zero between states derived from single group states 



^± 



tional to j (iSc^ — 9ciC2 + ISc^). It is not then possible to of Fj , T^ and Tf IRs as dictated by symmetry, selection 

rules permit their existence under double group formulation 



equate the difference in z components of angular momentum 
to integer multiples of h. 

In obtaining generators for irreducible perturbation trans- 
forming according to Fj symmetry, Luttinger used the two 
linearly independent matrices which form the matrix repre- 
sentation from Jp and Jf^ operators using a single set of 

basis^. The reason that J^ and JiJ operators return differ- 
ent linearly independent matrices, even when evaluated with 
respect to the same set of basis, is that they are different op- 
erators that have the same transformation properties. The two 
different linearly independent matrices obtained from a single 
set of adapted double group basis reflect the properties of two 
different operators, not the properties of the basis. Hence, the 
total angular momentum operator, is more appropriately de- 
scribed by Eq.dSll or Eq.(j9]l depending on the use of single or 
double group formulation. If the zone centre states can be de- 
scribed by adapted double group bases under the single group 
restriction, then the second invariant C2 for angular momen- 
tum operator would always be zero. The Zeeman interactions 
are different from the total angular momentum operator, and 



they are discussed in section VI 



The same issues arise when derving the matrix representa- 
tion of momentum operator p, or perturbation operator tt, be- 



tween states with Fg symmetry. As shown in Eq.( B5 1, it is de 



scribed by two linearly independent matrices under the double 
group selection rules, but only one matrix if adapted d ouble 
group bases are used under the single group restrictiorU^E^ 
Given the finite q Luttinger parameter observed in experiment, 
the use of adapted double group bases under the single group 
formulation is not appropriate. 

The cause of the incompleteness of the adapted double 
group bases under the single group restriction arises formally 
from the absence, or neglect, of the inter-band k independent 
spin orbit interaction. It is well known that spin orbit interac- 
tion will cause mixing between Fg IRs of the same paritjEISI. 
To properly consider the effects of perturbations to second or- 
der in fe, such mixing should be included. Then the zone cen- 
tre states of F^ symmetry in the conduction band would con- 



tain a mixture of bases derived from T^ 



F4 and F+ ® F3" 



IRs. The restriction to one single group IR is in principle at 
least, incorrect. One may argue that spin orbit effect is a small 
perturbation and such mixing would be minimal. However, 
additional mechanisms of mixing exist if one considers the 
way in which hybridised orbitals are constructed with con- 
sideration of spin (see section IVli. This mixing effect can be 



readily seen in the analysis shown in sections IX and XI 



The incompleteness of adapted double group bases when 
used under the single group formulation implies that the treat- 
ment of spin-orbit interaction is inappropriate from a more 
fundamental view point. Mixing may arise as a result of many 
electron effects/configuration in addition to the effects of the k 
independent interband spin orbit interaction. A consequence 
of this is that zone center states with F^ and F^^ symmetry 
may also experience the effects of spin orbit interaction as 
a shift in their energy, though the degeneracy of such states 
would not be lifted. While the matrix element of H^„^ is 



between states associated with F^, Fg and Fg IRs. The con- 
sequences of this observation are discussed in sectionMunder 
the double group formulation. 



IV. COMPARISON WITH 6-BAND HAMILTONIAN 
OBTAINED UNDER SINGLE GROUP FORMULATION 

The purpose of this section is to establish the relations be- 
tween the form of the 6-band Hamiltonian expressed in EWZ, 
and those in the literature obtained under a single group ap- 
proach using either perturbation theory or method of invari- 
ant. This should ease re-interpretation of any experimental 
data that relies on labelling of Mj numbers, its correlation to 
z component of angular momentum, and its use in designa- 
tion of heavy and light hole states. It should also highlight the 
difference in sign between the 72 Luttinger invariant in the 
Hamiltonian of EWZ, and those used in the literature. 

It is first important to clearly define the concept of heavy 
and light hole bands. The heavy hole (light hole) band is gen- 
erally interpreted as the set of Kramer degenerate pairs in the 
valence band that has the larger (smaller) effective mass at the 
zone centre in the (001) principle axes directions. The Mj 
number, serving as an index to the degenerate states at the 
zone centre, defines the order of the basis in the IR, and the 
form of the representation matrices of the Oh group. It is fre- 
quently correlated to z components of angular momentum. As 
stated in appendix |F] a direct link to angular momentum can 
not be made for point groups and any such link has to be based 
on the fact that Oh group is a subgroup of 0(3), and the bases 

of of and Df IRs of the 0(3) group also form the bases of 

2 2 

Fg and Fg IRs of the Oh group. Thus, the order of \J, Mj) 
bases of IRs of the Oh group must be the same as those in 
the \j,mj) bases of Dj^ IR of the 0(3) group, for the link 
between Mj and z component of angular momentum to ex- 
ist. For this to occur, the representation matrices of the 0(3) 
group differ in sign between proper and improper rotations. If 
the correlation between the AIj number to z component an- 
gular momentum is to be made, the bases of O/, group must 
be ordered in such as way that its representation matrices fol- 
lows the same rule. The double group bases used by EWZ and 
those given in appendixlAlfollows this convention. 

Within the 4-band model, the | J, Mj] 
scribe states with an effective mass along the (001) direc- 
tions given by '"^ , while the \J,Mj) = |f ,±5) pairs 

describe states with an effective mass of — ^^ — . The I J, MA 

71+272 I ' J' 

quantum numbers have the significance defined in section III 



I 2" , ± 5 ) pairs de- 



of EWZ, and determine the time reversal and spatial trans- 
formation properties of the specific basis function within the 
IR. Clearly, the Mj index of the heavy and light hole bands 
exchange with each other if the sign of 72 changes. Once 
the sign of 72 is determined by the material parameters (see 



Eq.( [D4b| and Eq.( [D5b] i in appendix|D]or Eq.(18b) of RefQ, 
the association of heavy or light hole bands with appropriate 



8 



Kramer pairs is determined. In appendix ID] a condensed ver- 
sion of the derivation of relations between the Luttinger invari- 
ants and second order interaction parameters is detailed, with 
a different choice of bases for the Fg IR as compared with 
those in EWZ. Since the definition of a second order inter- 
action parameter ensures that F, G, Hi, H2 are negative and 

C's are positive (except C+V+/^ '■^^ ^iS'^ °f 72 is mostly 
determined by the relative dominance between F and Hi in 
the single group formulation (see Eq.(D4bi), or (^^l _^ and 



^ s 1^ a 



C„l „+ in the double group formulation (see Eq.(D5b 1). The 

I 8 '1 8 I I 

sign of 72 is therefore mostly determined by the symmetry of 
the lowest zone centre conduction band states. In Ge and most 
other III-V compound semiconductors, a negative 72 may be 
expected due to the symmetry of lowest conduction band zone 
centre states being Fy (Fg), and the bases in F^ ordered in the 

same way as Dt- 

2 

The designation of valence states does not always corre- 



late with numerical value of experimentally measured effec- 
tive mass in different crystallographic directions. For exam- 
ple, when the in-plane mass is measured using cyclotron res- 
onance with the magnetic field applied in the [001] direction, 

the 4-band Hamiltonian yields a mass of — ^ — where the + 

J 71 ±72 

(— ) sign refers to the ||,±|) (||,±i)) pair The relation 
between these measured masses and Mj numbers would be 
in opposition to the heavy and light hole band assignments^) 
The energy eigenstates with an arbitrary k are also generally 
a mixture of zone centre states due to presence of off-diagonal 
terms in 4-, or higher-band models. In addition, one can not 
make an association between the mass measured in an arbi- 
trary orientation, with zone centre states defined with respect 
to the principal axes. 

The association of the heavy and light hole bands with the 
bases of the 6-band Hamiltonian, derivable from appendix D 
of EWZ, depends on the sign of Luttinger invariant 72. The 
Hamiltonian has been obtained from perturbation theory (see 
appendix [D]) and the method of invariant, and may be written 
as. 
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where. 



P - 71 fc' 



72(fcff-2fc2), 5 = 2%/373fc_fc, 



R = -Vii2{kl - kl) + 2^/?,i-i^k,ky, 



Rvs = -^iTikl 



2V3i 



k\\ = k^ + k , k — k^ 



P. 



llk\ 



If 72 > 0, the heavy hole states near the zone centre are as- 
sociated with \J,Mj) — ||,±|) basis functions. If 72 < 0, 
then the heavy hole states near the zone centre are associated 
with the ||, ±i) basis functions. 

Any unitary transformation can be made to the basis trans- 
forming as a F^ IR without affecting the dispersion relations 
of the valence bands. Thus, a unitary transformation can be 
considered which permutes the Mj = ±| and Mj = T5 
states. This is, in some way, equivalent to using the bases of 
Fg IR in place of the F^ IR. Most results in the literature were 
obtained using the ordering of the odd bases, even though the 
zone centre valence states of Ge must have positive parity, and 
those in most III-V compound semiconductors are dominated 
by the spatially even part of the wave functions. The group 



theoretical reason behind this unitary transformation is given 
in appendixlF] The transformation of the basis functions of the 
valence states are then describe by the relation, 

{Wi) ,W2) M'^) M'^)) ^ {\(^i) M2) M^) Ma))u 



c/ = 



/o -1 o\ 

1 

10 

\Q -1 0/ 



(12) 



where the basis transform as a row vector The U matrix has 
the same form as the T matrix originally identified in EWZ 
and corrected in Eq.(2) of the erratum^^. The T matrix trans- 
forms the basis of the Fg^ IR to ensure that the representation 
matrices for elements of the Fg IR of the Td group are ob- 
tained under compatibility relations. However, the U matrix 
is applied to the bases of F^ IR transforming as a row vector, 
so that U is actually the inverse of T defined by EWZ. The 
unitary transformation of the valance band and spin split off 
band is then given by. 



M 



U 04x2 
02x4 l2x2, 



(13) 
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The 6-band Hamiltonian referring to the transformed bases 



may be obtained from Eq.( 1 1 1 using the transformation in 
Eq.([T2]i and is written as. 
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-72. This leads to the inversion of the sign in front of the Q, R 



eff 



where terms have the same definition as in Eq.( 111. 

An effective Luttinger invariant can be defined as 73 
and R* terms, replaces R with R*, and R* with R. When a similar definition of effective Luttinger invariants, 72 

and 73 -^ — — 7I'* in the interband block is made, the sign on all terms in the interband block inverts. The form of the Lowdin 
contribution to the Hamiltonian is then the same as the one obtained by Luttinger and Kohri^, Suzuki and HenseP^I and Chao 
and Chuang221 and given by. 
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where, the definition of the terms differs from those in Eq.( 11|14| as, 

Q = j^/f{kl-2kl) = -j,{kl-2kl) 

R = -V3-t2^^{kl - kl) + 2VSi"f3k^ky = ^372(^:2 - k^) + 2VSi-f3k^ky 

Qvs = 72^^'(^| - '2kl) = -72^^(fcH - 2fc') 
Sy, = 2V5j3^^''k_k, = -2^37^ fc_fc^ 

-ttys 



Vil^f {kl-kl)+2Viill 



eff' 



fvxrvy 



VSlTikl - kl) - 2ViiiTkxky. 



(15) 



For the purpose of comparison, an approximation is made 
to restore the single group selection rules and ignore the spin 
orbit splitting in the definition of Lowdin terms. This approx- 
imation is consistent with those considered under the single 
group formulation '■*. Thus, utilising the single group F^ IR 
heritage of the 6-dimensional F^ F^ subspace, one can 
write the interband Luttinger invariants as. 



eff' 

72 = 

eff' 

73 = 



-72^ = 
-73^^ 



eff 

-72 =72 



73 



7i =7i 



(16a) 

(16b) 
(16c) 



This leads to equality among effective Luttinger invariants in 
the interband block, and their counter parts in the main va- 



lence band block in Eq.( 15 1. In this case, all the newly defined 
effective Luttinger invariants are positive, for most semicon- 
ductors with diamond and zincblende lattices, except diamond 
and possibly Si as indicated by Lawaet^^. 



The form of Hamiltonian in Eq.( 15 1 is similar to those ob 



36 . 



tained by DresseUiauP, Kans^i and Luttinger and Kohn 

ing the either the alternative diagonalisation matrix or bases 

functions of the Fg and F^^ IRs. Any difference in the signs 



of the terms between expressions in the literature and Eq.( 15 1 



is due to different choice of phase in the bases. The form of the 



Hamiltonians shown in Eq.( 1 1 1 and Eq.( 15 1 are related to each 



other by the unitary transformation given by Eq.( 13 1, and are 
therefore both correct as they give rise to the same dispersion 
relation, for a given set of Luttinger invariants. However, the 
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the methodology emp loyed to obtain Eq.( 15 1 using the basis 
of {Tq, rg_jrg }li2136lcan only be justified under the method 
of invarianP^l as shown in appendix^ The basis of Hamilto- 
nian shown in Eq.( 15 i, remain as those of the F^ffiFy IRs. As 
shown in appendix IFI the association of Mj with the z com- 
ponent of angular momentum is not correct for the form of 



Hamiltonian in Eq.( 15 i, and the subsequent association with 
heavy and light hole band (accepting z components of angular 
momentum as Mjh) is not valid. 

The zone centre energy eigenstates in the crystal are the 
eigenstates of the angular momentum operator but the eigen- 
values for J^ and j^ are not generally integer multiples of 
fi^/A and h/2 respectively. Any link between the Mj quantum 
number and z component is, at best, a correlation. In dealing 
with physical effects such as optically induced transitions, it 
is important to evaluate the relevant matrix elements in the 
Fermi Golden rule that determines transition probability. For 
example, the relevant perturbation operators, for transition in- 
duced by light polarised as cr± in the Faraday configuration, 
are p^ ± jpj^ll The permitted transition can easily be ob- 
tained using the form of M^„ ^^ ± iM^„ ^„ matrices given 
in appendix Ib] For light linearly polarised in the (001) direc- 
tion, the corresponding transition matrix between states of the 
valence band (F^), and lowest conduction band (F^ ) is given 

by, 



M 



p+ p- 
^ 8 '^ 7 



_?z. /O 2^ 
m I 2 Oy 



It is clear that only states ||,±|) in the valence band are 
coupled to the conduction band states with Fy symmetry by 
light polarised in the (001) direction. As pointed out in the 
erratum to EWZ and earlier, the nature of | J, Mj ) — || , ± |) 
(heavy or light) depends on the sign of 72. For materials with 
negative 72, these states corresponds to the light hole band as 
observed in Ge and many other compound semiconductors. 



The unitary transformation indicated in Eq.( 12 1 must also 
be applied to other interaction matrices, such as the strain de- 
pendent perturbation, and as a result the sign of the defor- 
mation potentials should be carefully considered. However, 
an important consequence of this impUcit unitary transfor- 
mation is the incorrect association of coupling provided by 
Qys term between the ||, ±i) states in the valence band and 
1 5 ' "F I) spin split off states, based on the form of Hamiltonian 



in Eq.( 15 1. Such an association has been used in the interpre- 
tation of experimental data of materials under uniaxial stress, 
and underpins some of our understanding of the the electronic 
structure in these crystals'**^ "*■'. 

In the rest of the manuscript, any reference to the 6-band 
Hamiltonian refers to Eq.( 111. Reference to Eq.( T4p5 1 are 
stated explicitly. 



V. DOUBLE GROUP FORMULATION 

It appears that all the three elements in the treatment of 
spin-orbit interaction under the single group formulation had 



some deficiencies. The use of double group formulation is 
then a natural way of dealing with spin-orbit interaction, by 
treating the k independent spin orbit term in Eq.(fT]l directly 
in the un-perturbed Hamiltonian. The discussions in section 
III indicated that the adapted double group basis were incom- 
plete when utiUsed under a single group approach. This was 
attributed to a restriction placed upon such bases by the lack 
of mixing under the single group formulation. In this section, 
mixing is considered under both the single and double group 
formulations to establish the key difference between the two 
approaches. Other impacts of the double group formulation 
on the analysis of material parameters is also considered. 

One of the main distinctions between double group for- 
mulation of EWZ and traditional single group formulations, 
is the use of hybridised orbital-derived zone centre states as 
bases for the expansion of Hamiltonian H{k), instead of those 
derived from atomic orbitals. The action of crystal field on 
the atomic orbitals leads to zone centre states with symme- 
try given in Table III] In constructing zone centre states from 
atomic orbitals, a degree of mixing of atomic orbitals is per- 
mitted under the single group formulation, provided that they 
are compatible under the single groufP^. The compatibility 
between atomic orbitals and single group IR are shown in Ta- 



ble III We will consider these mixing with the condition that 
they do not alter the angular dependence of first order fe • tt in- 
teraction under the single group formulation. In other words, 
mixing of atomic orbitals is subject to compatibility of sin- 
gle group IRs. Thus, states with F^ Fy" symmetry derived 
from the single group F^ IR contain contributions from p, dxy 
and /^3 orbitals, but no contribution from ^^2 or fz{x^~y^) or- 
bitals. Under the single group formulation, incorporation of 
zone centre states with F^ (1^3^) and F^ V^ (F^) symme- 
try containing d^^ and fz(x^^y^) orbitals, respectively, would 
lead to new bonding states which are not generally observed 
in experiment. They therefore must be considered as non- 
participating in the formation of the covalent bond. Likewise, 
under the single group formulation participation of f*ix2^y2\ 

orbital in anti-bonding states with F^ symmetry is not pos- 
sible, as this orbital is derived from a single group IR with 
F^ symmetry instead of Fj" symmetry. It is clear that un- 
der the single group formulation, restrictions are placed upon 
the mixing of different symmetry compliant orbitals into the 
bonding and anti-bonding zone centre states. 

Under the double group formulation, the action of crys- 
tal field and spin orbit interaction is considered together as 
part of the un-perturbed Hamiltonian. Full mixing of orbitals 
within individual double group IRs is permitted. This can 
be explained under the effect of hybridisation under the ac- 
tion of both the crystal field and spin orbit interaction, as in- 
dicated in Table 



IV 



Bonding states with Fg symmetry and 
anti-bonding state with F^ symmetry, may then arise from 
p,dxy,dz2fz3,fz{x^^y^) and s* , f*^^, f*^^2_y2-j orbitals, re- 
spectively. The fc-independent spin orbit interaction leads 
to further mixing between some states with the same parity 
in crystals with the diamond lattice, and among all states in 
crystals with zincblende lattice. All the atomic orbitals are 
permitted to participate in the formation of bonding and anti- 
bonding zone centre states. In other words, the double group 
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TABLE II. Atomic orbitals symmetry and their contribution to bonding and anti-bonding states 



Atomic 
orbital 

s 

P 



Symmetry 



^ 2 (jxyz) 

r4(/.3) 



Bonding states 



No spin 






With spin 



^ 8 



Antibonding states 



No spin 



With spin 



r^(-) 


r^ 


r^ 


r^- 


Tf 


^4{P) 


^t 


r+®r+ 


r4" 


^S ® ^6 


|r+(d.,) 


rj 
r^ 


r+ffir+ 

r+ 

-•- 8 


r4 


Fg- ® Fg- 



J7 

Fg ® Fg 

-L « © I 7 



TABLE III. Mixing of atomic orbitals permitted under single group 
formulation 



Symmetry (zone 
centre state) 



Contributing atomic orbitals 



^^ 
^^ 

Fr 



.r+ 



'r6+ 



s, 7x1/2 
s* f* 

"^ ■) Jxyz 
J z{x^ ~y^) 



TABLE IV. Mixing of atomic orbitals permitted under double group 
formulation 



Symmetry (zone 
centre state) 



Contributing atomic orbitals 



f; 



S, Jxyz, Jz{x^—y^) 

Pi dxy, J z^ 

P, Clxy, Clz^, Jz^, Jz(x^—y'^) 

P*, dly, f*3 

P*, dly, d*2, /*3, /*{^2_j,2) 

S ) Jxyz, .fz(x'^-y^) 



formulation calls for symmetrization of linear combination of 
wave functions to be performed at double group level, and the 
Hamiltonian of the unperturbed system includes the k inde- 
pendent spin orbit interaction. 

The energetic ordering of the zone centre states has a strong 
impact on the sign of the 72 Luttinger invariant in the valence 



band, as stated in section |IV] and shown in appendix |D| This 
ordering also has a strong impact upon the bandstructure of 
the semiconducting material. From the point of view of a 



semi-empirical technique, the energies of zone centre states 
are generally determined from experimental observation. As 
a result, mixing allowed under the double group formulation 
would not alter the the assignment of zone centre state ener- 
gies in practical implementation of the k ■ p method as they 
are determined from experiment. However, under the double 
group treatment of spin-orbit interaction, the effects of spin 
can shift the energy levels of states with Fy and F^^ symme- 
try. This provides an additional mechanism through which 
the ordering of zone centre states in the conduction band can 
be modified. Without solving the unperturbed Schrodinger 
equation under the double group formulation, group theo- 
retical methods indicate that spin can effect states with F^ 
and Tq symmetry in addition to the lifting of degeneracy be- 
tween states with F^ Fy" and Fg" Fg" symmetries. Since 
WV{r) X p and S both transform like an axial vector (T^) 
under the Oh group, the H§^i term transforms according to 

) rt, and the interaction is time 



.F+. 



>^^ 



reversal even. The product decomposition of states with either 
F^ or V^ symmetry both give [Fj'"] {F^} IR in the decom- 
posed direct sum. Thus, selection rules indicate that states Fy 
and Tg symmetry (and states with F^, Fg", Fg ) contain at 
least a [T^] IR in their decomposed direct sum from their di- 
rect product, permitting spin orbit interaction to shift the en- 
ergy levels of such states at the zone centre. Consequently, 
many electron effects should be included in the determination 
of zone centre state energies or treated as a fitting parameter 
under the semi-empirical scheme. 

It is particularly interesting to examine the trend from di- 
amond to a-Sn, which indicates a reduction in energy of the 
F^ state relative to other states. In materials with small spin 
splitting effects (diamond and Si) the Fy state resides at en- 
ergies above the Fg^ Fg^ states. In Ge, with large spin split- 
ting, the Fy state has shifted below the Fg^ Fg states and 
in a-Sn with the largest spin splitting inversion has occurred. 
The double group formulation may then provide an explana- 
tion why states with F^ symmetry reside at lower energies 
than states with Fg © Fg symmetry in materials such as Ge 
or GaAs but not in materials such as diamond and Si. Without 
spin-orbit interaction, the lowest F^ basis has I = 3, whereas 
the lowest Fg © Fg bases have 1 = 1. While this may ex- 
plain why the Fy state should have higher energy among the 
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degenerate anti-bonding hybridised orbitals (higher number 
of nodes associated with larger I number indicate higher ki- 
netic energy), as in the case of diamond and Si, spin orbit 
interaction should be taken into account to decide the final 
ordering of the conduction states. With the large orbital an- 
gular momentum quantum number associated with the bases 
of r^ IR, one could expect a larger downward shift of F^ 
states due to spin orbit interaction compared with that of the 
Fg Tq states. This is, of course, dependent upon the devia- 
tion of zone centre states with F^ symmetry away from those 
prescribe by the single group formulation (participation from 
s* , /* 2 orbitals only, as shown in Table llli. The mixing, or 
introduction of f*,„2_, 2-, orbitals into states with Fy symme- 
try under the double group formulation, provides the required 
deviation from a single group prescription. Thus, in the case 
of semiconductors comprising of heavier elements, such as Ge 
and a-Sn, one can therefore expect a stronger spin orbit effect 
and a larger degree of mixing in the zone centre states. Spin- 
orbit interaction and mixing may then have a deciding effect 
on the ordering of the conduction band zone centre states in 
these materials, even though no such effect is expected under 
the single group formulation of the k • p method"*^ . 

The effects of mixing in zone centre state as described un- 
der the double group formulation may be illuminated further 
by making a comparison between crystals with a diamond and 
zincblende lattice. In crystals such as SiC or AlP, the conduc- 
tion band state with Fg symmetry is below states with Fg and 
F7 symmetry. By considering the compatibility relations be- 
tween the Td and Oh groups, these materials behave like Ge 
or GaAs, in terms of the ordering of zone centre states in the 
conduction band, even though the constituent elements belong 
to the same row as, or above that of Si in the periodic table. 
The distinction between them is the lack of inversion sym- 
metry in crystals with zincblende lattice. This allows the k 
independent spin orbit interaction to couple the valence band 
zone centre states to those of the conduction band under the 
double group selection rules. This additional coupling, over 
what was allowed in crystals with diamond lattice, enhances 
the deviation of the states in the conduction band with Fg sym- 
metry from those prescribed by the single group formulation. 
Therefore, spin orbit interaction and the additional coupling 
expected in zincblende crystals such as SiC and AlP, can en- 
hance the downward shift of states with Fg symmetry towards 
the valence band in comparison to those with F7 and Fg sym- 
metry. 

In the latter part of this section we shall examine the impact 
of double group formulation on the analysis and interpreta- 
tion of material parameters. The sign of 72 for Sp^ is gener- 
ally considered negative within the double group framework 
of EWZ, as in the case of Ge, GaAs and many other semi- 
conductors. To determine if this is the case for Si, one could 
measure the ratio of heavy and light hole masses in SiGe al- 
loys. If the Luttiger invariants in Ge and Si have the same 
sign one could expect the heavy/light hole mass ratio to be 
a monotonic function of composition from Ge to Si, if the 
Luttinger invariants are smooth monotonic functions of com- 
position. This is indeed the conclusion of simulation based 
on Kane's model reported by Braunstein'*^. However, the ex- 



perimental data shown in Fig. 10 of Refl46l indicates that the 
mass ratio approaches unity, but would then clearly have to 
increase back to the value of Si. This indicate a point of in- 
flection exists in the experimental data. It can only occur if 
the value of 72 for Si is positive and the point of inflection 
corresponds to 72 = at the given composition. Therefore 
the experimental data here suggests that Si, like diamond, has 
a positive 72, while Ge has negative 72. In the SiGe material 
system, one must therefore deal with possibility of 72 chang- 
ing sign depending on composition. The scaling of various 
material parameters under envelope function approximation 
should take this into account. 

The question is then raised about the interpretation of data 
by Hensel and Feher^, and Hasegawa"*'. Taking into account 
the difference in the order of basis as detailed in section |IV] 
the determination of sign of 72 relies on application of specific 
second order perturbative term which leads to linear stress 
dependence on effective mass of hole states coupled to the 
spin orbit band. However, it is easy to see that the effective 
mass of all states in the valence band would be dependent on 
the applied stress (products of strain and quadratic k terms 
in the Hamiltonian contains F^'" irreducible component) from 
method of invariant. This term has the same order as those 
considered by Hasegawa^. The conclusion that stress affects 
the effective mass of only one of the valence band states may 
not be correct. One aspect of the measurement of Hensel and 
Feher shows very similar behaviour of Si compared with Ge 
and most other semiconductors with diamond or zincblende 
lattice. Under uniaxial compression in the (001) direction, 
the states with lighter in-plane mass (normally referred to as 
heavy hole band) are depopulated. Thus, the energy of the 
heavy hole states decreased relative to the energy of the light 
hole states under uniaxial compression. This would be con- 
sistent with Si having a negative 72 as in the case of Ge and 
GaAs, if the deformation potential b has the same sign in both 
these materials. This result would then be in contradiction to 
results of Braunstein"*^. 

Similar optical reflectance measurements under uniaxial 
stress exist in the literature and should be able to resolve this 
issue, at least in principle. The experimental data for GaAs 
and Ge"^- "*^ are relatively easy to interpret. There is clear evi- 
dence that the light hole states move upwards in energy under 
the uniaxial compression in the (001) direction and its value 
has a non-linear dependence on the stress. This leads to the 
conclusion that 72 is negative and heavy hole states are as- 
sociated with \J,Mj) = ||,±|) quantum numbers. Data 
for Si is inconclusive42^. The quadratic dependence asso- 
ciated with \J,Mj) = I I, ± 1) states can not be identified 
clearly from Fig.4 of Ref 43 This may be due to the small 
spin orbit splitting in the valence band of Si compared with 
the large strain perturbation actually applied over the range of 
stress used. Judging from data of Hensel et al'*", decoupling of 
states occurs at stress of approximately 4 x 10^ Pascal. This is 
smaller than the smallest none zero stress of 10^ Pascal used 
by Laude et al.™ Therefore, the sign of 72 for Si remains un- 
determined. 

The relations between Luttinger invariants and second order 
interaction parameters are important when obtaining material 
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parameters for the k • p method. Thus it is important to re- 
late all the Luttinger invariants to the second order interaction 
parameters, under both the single and double group selection 
rules for appropriate comparisons of the two formulations. In 
absence of magnetic field, there are only three Luttinger in- 
variants 7i, 72 and 73. To extract the five second order C, pa- 
rameters defined in EWZ from the experimentally measured 
Luttiger parameters, would then be an under determined prob- 
lem. With magnetic field applied, two further invariants k and 
q are available in addition to 71 , 72 and 73, ensuring that there 
is a one to one correspondence between invariants and second 
order interaction parameters. Under the DKK model, there 
are a total of four invariants, L, M, N and the magnetic in- 
variant P, which can be related to a set of four second order 
interaction paramters, F, G, Hi and H2- Efforts to relate 
the five invariants obtained under the double group selection 
rules with second order interaction parameters in other models 
based on single group formulation lead to linear dependence 
between the invariants'^ ■'^■'^"*'*, contrary to the symmetry ar- 
guments. In the following section, a full set of relations are 



extracted between invariants and second order interaction pa- 
rameters under the double group selection rules by consider- 
ing the magnetic perturbation. 



VI. ANTI-SYMMETRIC LOWDIN TERM 

In this section the anti symmetric Lowdin term is derived 
from perturbation theory and group theoretical method under 
the double group formulation. The relations between all the 
Luttinger invariants and second order interaction parameters 
are derived. This enables a one to one mapping between the 
Luttinger invariants and independent second order interaction 
parameters in each block of the Hamiltonian. Representation 
of Zeeman interactions are then obtained from relevant terms 
in Eq.([T]i. 

The matrix representation of Hamiltonian between zone 
centre states associated with representation a" and a" (clas- 
sified according to double group representations) in the pres- 
ence of magnetic field may be written in block form as. 



_fiQU qu — _rzQ 



a ,Q , rrCt ,Q 






<y-E^/^^""." 



(17a) 
(17b) 
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x,y,z x,y,z U 
Ldwdin:S / j / j / j 



rjOL ,a 

Ldwdin:A 






ky{Ec^ + Ec^ — 2Ec,iii) 



x,y,z x,y,z U 

EEE- 

Hso3 — gol^sSa^^a-" ' B 



K{E^.+E^.~2E^f<) 






4(i?a" - E^f3){Ea^' - E^ 



(17c) 

(17d) 
(17e) 



where Af^u ^i. are defined in appendix [bI H^ '" is diago- 
nal and contains zone centre state energies when a" = a", 
H? '" is the first order k ■ n interaction, and i?? '", is the 
interaction between states belonging to a" , a" IRs mediated 
by remote states in a^ IRP, and is partitioned into symmet- 
ric and anti-symmetric parts. The Zeeman interaction due to 
electron spin (Hsos) arises directly from the Foldy-Wouthysen 
transform. The anti-symmetric Lowdin term arises from lack 
of commutability between components of wave vectors due to 
presence of magnetic field. The total Zeeman interaction con- 
sists of contributions from the anti-symmetric Lowdin terms 
and those due to electron spin (Hgos)- 

With second order interaction parameters as defined by 
Eq.(l 1-13) in EWZ, the symmetric part of the Lowdin term is 



obtained from Eq.( 17c 1 using M matrices shown in appendix 
[B] The form of this symmetric Lowdin term has been given in 



Eq.( 111, and relations between Luttinger invariants and sec- 



ond order interaction parameters are given in Eq.(D5 1. 

For the anti-symmetric part of the Lowdin term, we first 
consider interactions mediated by the remote states with Tg 



and Fy symmetry. Constructing the anti-symmetric Lowdin 
term from Eq.(17di using M matrices obtained in appendix 
iBlyields two distinct terms. 



Lowdin-.A 



ih 






H 



Ldwdin:A 



c.p.)(18a) 

-86 

-2AiBCpT p+ - {Jx[ky, fc.] + c.p.)(18b) 



^P+ r 

1 a ,1 



8 e 



where c.p. stands for cyclic permutation of the preceding term, 
[fc^, k^] = k^kv — k^k^, and the two generator matrices J and 
J are the same as those used in Eq.( |7|9[ ). The form of these 
matrices and additional phase factor "i", are obtained from 
constructing the anti-symmetric Lowdin terms using pertur- 
bation theory. Both sets of generator matrices transform ac- 
cording to rj IR of the Oh group and are clearly linearly in- 
dependent of each other. As in the case of symmetric Lowdin 
term shown in EWZ, the interaction between states with Fg 
symmetry contribute three distinct interaction parameters to 
the anti-symmetric Lowdin term. 
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H 
H 
H 



Ldwdin:A 

r, 



F~ 1 tri v~ 1 111 

^J.BCl p+ — {Jx[ky, k^] + c.p.) + /isCp" p+ — 



yJx[Ky^ kz 



\in:A = 4Mi3Cpt p+ — (Jx[fcy, fcz] + C.p.) 



Ldwdin:A 
L6wdin:A 



TJ.n 



c.p.) 



= 5Mi3Cp+'p+— (Jx[fcy,fcz] +C.p.) +3AiBCp+ p+— ('/a^ifcyjfcz] + C.p.) . 



(19a) 
(19b) 
(19c) 



r 



Application of a magnetic field in any of the principal 
axes direction leads to the consideration of the problem un- 
der envelope function theory as prescribed by Luttinger and 
Kohn^^. The coupled differential equation is obtained by re- 
placing only the components of wave vector in the Hamilto- 
nian perpendicular to the magnetic field with fc„ — > A;^ = 
jg^ — h f^a, where e is the positive universal constant of 
electronic charge. For magnetic field aligned in the z direc- 
tion, the commutation relation between the corresponding op- 
erators is, 



[kx^ ky\ 



eB 

ih 



(20) 



under the symmetric gauge (A^ = — 2^ , Aj, = ^ , yl^ = 0). 
The wave vector component k^ remains a good quantum num- 
ber as translational symmetry along the direction of magnetic 
field is maintained. 

The matrix representation of Hgos with respect to the single 
group product bases is well defined since the spin and orbital 
parts are separated. A unitary transformation into the adapted 
double group bases yield a contribution to the Zeeman inter- 
action with modification to the k parameter by a factor —go/6 
as shown in Eq.( D4 1. The problem is more complicated given 
the mixing of the zone centre state. However, S transforms 
according to T^ IR of the O;, group. Therefore, the matrix 
representation of S with respect to the double group bases of 
r^ can also be expressed in a similar way to Eq.J9Jl, but with 
different set of reduced tensor elements. 



■J^ — "■ (l^aOfj, + QaJfj, 



(21) 



where k^. and q^^ reflects the properties of zone centre states 
with r^ bases due to mixing. As generators of transforma- 
tion matrices for spinor states under rotation operation, S'p 
are subject to the same constraint as the angular momentum 
operator shown in Eq.( 10 1. It is not clear how k„ and q^ cab 
be determined from theoretical or experimental perspective. 
Writing the Zeeman interaction as. 



^Zeeman — ^ Ldwdin:A ^~ ^so3 

= -2fj.B [k {JxBx + c.p. 



and combining results from the symmetric part of Lowdin 
term, the associated Luttinger invariants are then given by. 
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The anti symmetric Lowdin terms can be obtained for the 
conduction and spin-orbit split off bands, as well as the in- 
terband block between the valence and spin-orbit split off 
bands using the same methodology. For conduction"*'' and 
spin-split off bands, the Zeeman interaction is written as, 

Hzeeman — HLowdin.A + H sod, = —l^f^sf^K {OxBx + C.p.). 

The Luttinger invariants for these blocks are related to cor- 
responding double group second order interaction parameters 

by, 
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(23h) 
(23i) 



q {JxBx + c.p.(^2) 



For the interband block between the valence and spin-split off 
band, the Zeeman interaction may be written as, H Lowdin, A ~ 
— /i^/iK^'' {Jx^Bx + c.p.), where J'"^ is written as, 
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The form of these matrices were obtained from perturbation 



theory using Eq.( 17di. They may also be obtained using the 
method of invariants by selecting the appropriate generating 
operator (irreducible tensor operator) along with the appropri- 
ate bases of T^ and Fy IRs of the O;, group, as detailed in 
EWZ. Alternatively, one may consider the total angular mo- 
mentum operator J evaluated with respect to the j ~ 5/2 
bases of F^ and F^^ IRs of the Oh group. The invariants are 
then related to the second order interaction parameter by. 



72 = 



73^^ 



^/2 (^5 
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^r+ r+ 
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(23k) 

- 30^^^(231) 



Thus, there is a correspondence between the Luttinger invari- 
ants and second order interaction parameters in each block of 
the Hamiltonian. These relations may be compared with those 
of the DKK Hamiltonian with the k independent relativistic 
correction detailed in section Vlll as shown in Eq.( |D4| i, or 
Table I of Ref 50 Under the single group formulation, the 
following relations are obtained between Luttinger invariants 
in different blocks of the Hamiltonian, 
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1^" = 50 




(24a) 
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With the adapted double group bases under the single group 
formulation, the contribution from Zeeman interaction due to 
electron spin is a constant matrix with k,„ = ^] q^ = 0, 
K% — — i, and k'"J' = ^. Thus the contribution from Zeeman 
interaction due to electron spin makes the Luttinger invariants 
associated with magnetic interaction different among the dif- 
ferent blocks even under the single group formulation. 

It should be noted that there is freedom in choosing the form 



and the scaling of the generator matrices shown in Eq.( 18 1 



Thus, the specific relations between Luttinger invariants and 
second order interaction parameters are dependent on the 
choice of linearly independent generator matrices. The choice 
of generators made here is guided by perturbation theory, and 
naturally leads to simpler relations between the Luttinger in- 
variants and second order interaction parameters. In the liter- 
ature, the Jfl generator is generally used as a second linearly 
independent generator and is associated with the q invariant. 
Considering the double group bases used in this manuscript, 
J^ is equivalent to J7^ = ^J^f^ + jJp,- This means that the k 
and q invariants defined in this manuscript are related to those 
generally quoted in the literature {k' , q') by the following re- 
lation. 



2^ 



3 , 



(25a) 
(25b) 



The treatment of magnetic interaction for B along the [001] 
direction yields results that are no different from those of 
Luttinger**. The distinction between J^ and Jf^, can easily be 
explained in terms of linear combination of generators. How- 
ever, there is a difference in the process of deriving Luttinger 
invariants from experimental data. If the magnetic field is ap- 
plied the the [001] direction, the g factors for the valence band 
are given by. 



511 (M, 
511 (M, 



-2k + 6(7 
-6k + 2(7 



(26a) 
(26b) 



These relations may be compared with the expression given in 
Table I of Refl48las. g\\ ^ 2k' + \ = 2k - 6q. In the Hter- 
ature, the magnetic invariants k' and q' are related to second 
order interaction parameters obtained under single group for- 
mulation using single group selection rules^^. The measured 
parameters are thus linearly related to one another, contraven- 
ing double group selection rules and introducing errors in their 
values. The mixing of zone centre states also requires a differ- 
ent treatment of contribution to Zeeman interaction from the 
electron spin compared to single group formulation. 

A point to note is the difference between matrix representa- 
tions of total angular momentum shown in Eq.(|9]l and Zeeman 
interaction shown in Eq.(22 1. While the two expressions share 



the same generator matrices, the reduced tensor elements are 
not necessarily related. The coefficient ci and C2 for angular 
momentum are subject to constraints discussed in section III 
whereas k and q of Zeeman interaction are not. Hence, the 
general relation between Zeeman interaction and total angular 
momentum is lost in the crystal. This distinction arises from 
the dependence of anti-symmetric Lowdin interaction on both 
the remote states and near states in contrast to the sole de- 
pendence on near state for angular momentum. This point 
is demonstrated by the results of anti-symmetric Lowdin term 
and angular momentum obtained from bases of F^ IR derived 
from r^ IR. There is only one generator for J, but two for 
the anti-symmetric Lowdin term evaluated from perturbation 
theory. As a consequence, the unitary transformation which 
diagonalise the Zeeman interaction in a given orientation de- 
fined by external field differs from that generated by the total 
angular momentum operator J through rotation of axes. 

In other orientations, the perturbation approach illuminates 
the origin of the anisotropic nature of the interaction associ- 
ated with q invariant, and allows the derivation of an analytic 
expression for its matrix representation for the valence band 
for arbitrary orientations. Magnetic interactions with the B- 
field directed in other orientations is discussed in the follow- 
ing section. 



VII. MAGNETIC INTERACTIONS WITH B-FIELD IN 
OTHER ORIENTATIONS 

With magnetic field applied in directions other than the 
(001) directions, appropriate changes must be made to con- 
struct the relevant Hamiltonian. In the context of method 
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of invariant, the p roced ure has been described by Luttinger^ 
and is used wideljHJH) ^jj alternative method based on the 
transformation of wave vector only, was used by Eppenga and 
Schuurmans^"*. The method of Eppenga and Schuurmans dif- 
fers from others in terms of the bases used. A perturbation 
theory based approach has been applied to study dispersion of 
crystals under uniaxial stress^ by projecting the bases states 
of {\X) , \Y) ,\Z)} in appropriate directions. An operator or- 
dered Hamiltonian for treatment of spatially confined systems 
was also obtained using the same approach^'* ^^. In this sec- 
tion, we derive the necessary unitary transformation induced 
by a change of the coordinate axes, and derive the Lowdin 
interaction from the perturbation theory approach. The differ- 
ence between methods of Luttinger and Eppenga and Schuur- 
mans are clearly identified. 

It is important to note that there are two coordinate sys- 
tems, and the associated orthonormal bases of vector space, 
spanned by the degenerate zone centre energy eigenstates, 
form the representation the symmetry group. The {x,y,z} 
axes are oriented along the (001) crystallographic directions 
whereby the z axis is chosen as the quantisation direction of 
the third component of angular momentum. The point group 
symmetry of the crystal defines this coordinate system. The 
orthonormal bases are simply the unperturbed zone centre en- 
ergy eigenstates referring to the {x, y, z} coordinate system 
I't'a^^i)- The second coordinate system with {123} axes, is 
defined by the externally applied magnetic field B, which is 
parallel to the 3 direction and defines the quantisation direc- 
tion of the third component of angular momentum. The cor- 
responding zone centre energy eigenstates refer to the {123} 
coordinate system \(t>^^i)- In principle, the Zeeman interac- 
tion should be diagonal with respect to such bases. In apply- 
ing the envelope function theory, as prescribed by Luttinger 
and Kohn^^, the replacement of wave vector components with 
operator occurs in the {123} coordinate system. Using the 
symmetric gauge, we have. 



{123} coordinate system in terms of wave vector components 
and optionally associated bases. The Lowdin interaction de- 



fined in Eq.( 17 1 contains a typical term. 



fc™A^". ^,Ml, „.K 



where m and n index the directions in the {123} coordinate 
system, and A^™, ^^ are evaluated with respect to the 10;!^??) 
bases. This typical term may be related to M"\ ^ referring 
to the \4u^J') bases, through the unitary transformation defined 



inEq.(28i 



fc™X" ^,Ml, „.fc„ = A:„,C/"" W™ „,3C/"'[/"' W", „„C/""fc„ 



[/" 'fc„M^„„,M^,„.A;„C/" 



Thus, one can evaluate Af' ' 



with respect to the 



jxyz^^ 



bases, then construct the symmetric and anti-symmetric 
Lowdin terms, and finally make the unitary transformation 
into the bases of \(b]?'?\. The interaction matrix, M™ « can 
be calculated using a projection of the momentum operator 
pi,P2,P3 along the direction {a;,y, z}. Thus, if the compo- 
nents of a spatial vector in the second coordinate system are 
related by. 



R\\ Ryi i?i3 
-R2I ^22 ^23 
i?31 i?32 i?33 




(29) 



where the matrix R is real and orthogonal, then the relation 
between interaction matrices is given by. 
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A typical term in the Lowdin interaction is then given by. 



fci — > fci = 



fc2 = 



fcl,fc2 



l_d_ 
i dri 

i_a_ 

i dr2 

eB 

ih 



eBr2 
2h 

eBri 
2h 



(27a) 
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The third wave vector component k^ remains as a good quan- 
tum number. It is not possible to define the corresponding 
commutator in the {xyz} coordinate system, since the com- 
mutation relation defined by Eq.(27c 1 only holds in the {123} 
coordinate system. The bases of the vector space in the two 
coordinate system are related by a similarity transformation. 
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k h 



Therefore a general term in the Lowdin interaction transforms 
according to the direct product of two I — 1 representations 
of the 50(3) group, which is the symmetry group of proper 
rotations in 3-dimensions. 

Let us now focus on the magnetic interaction and anti- 
symmetric Lowdin term. In applying these relations, we ob- 
tain exactly the same results as in Eq.( p"8|23 i with the wave 
vector component referring to the {123} coordinate system, 
but the corresponding generator matrices defined still refer- 
ring to the bases |0qX^). Thus, the Zeeman interaction can be 
written as. 



Ceenian 



'2fiB {nJ^^ 



qJDB 



(32) 



In principle, the action of this unitary transformation should 
diagonalise the Zeeman interaction. 

Prior to the procedure of formulating the envelope function 
theory, and subsequent second quantisation of Landau levels, 
one must obtain the un-perturbed Hamiltonian referring to the 



The invariants, as one may expect, are un-affected by the co- 
ordinate transformation, and their relations to the interaction 
parameters remain as described by Eq.(23i. The generator 
matrices, on the other hand, also transform according to the 
I — 1 representation of the 50(3) group. This arises from 



17 



the fact that the anti-symmetric part of the product representa- 
tion also corresponds to the I — 1 representation of the 5*0(3) 
group. For the valence band, we then obtain two sets of lin- 
early independent generators, 

•^m'^ = RmlJ'x + Rml'Jy + RmsJz (33a) 

J^m^ = RmlJx + Rm2Jy + Rm?,Jz (33b) 

where the superscript xyz indicates that these generators are 
referring to the basis |0q«^), and the unitary transformation is 
still to be performed. They can be obtained either from these 
direct projections using Eq.( 



where J^ always has the diagonal form, but J3 is defined by 
Eq. p5a[ ). The J3 generator is equivalent to the 6*3 matrix in 
Ref 8 , and agrees with the form used in the paper (see Eq.(58)) 
once the differences between generators is taken into account. 
For the Fg IR, the corresponding matrices can be obtained 
by using the unitary transformation defined by the matrix T in 
Eq.(G3 1. In other words J^'^'^ will be diagonal and J^'^^ will 



be anisotropic. In the other interaction blocks of the Hamil- 
tonian, the Zeeman terms are diagonalised into the same form 
as in the {xyz} coordinate system. 



We note that the unitary transformations defined in Eq.( 28 1 



Eq.( 17di with aid of Eq.(31 1. The only relevant terms in the 
anti-symmetric Lowdin term are those related to [fci, k2\, as 
the magnetic field is along the 3 axis. The other two terms 
are zero since ^3 remains a well defined quantum number and 
commutes with ki and ^2- 

To re-quantize the zone centre states along the 3 axis, one 
needs to diagonalise the Zeeman interaction, shown in Eq.(|32]i 
for the valence band. In general, the two linearly indepen- 
dent matrices can not be diagonalised simultaneously, except 
when the 3 axis is in the (001) directions. The required 
unitary transformation of bases, given in Eq.([28]l would be- 
come dependent on the material parameters k and q and leave 



33]), or constructed directly from or those used in its place and defined by Eq.(|34]l are not neces 



xyz\ 



7123 



123 



qJz 



123 



) as a whole diagonal, but both J^ and 
J^^"^ individually non-diagonal. It is inconvenient and a uni- 
tary transformation which diagonalise J is used. It was mis- 
takenly suggested that this alternative transformation quan- 
tises the degenerate states split by Zeeman interaction in the 
3 axis*^, but this is not true since J3 remains non-diagonal. 
The unitary transformation which diagonalises the J genera- 



tor and is used in place of Eq.( 28 1 is given by, 



u^i 



cxp[-ij' • 6] 



(34) 



where 6 defines the rotation of coordinate systems from 
{xyz} to {123}. The transformation matrices for magnetic 
field aligned in [110] and [111] directions have been obtained 
in appendixO The diagonalised J^^^ generator has the same 
form as J^. The form of the Jl^^ generator for the magnetic 
field aligned in [110] and [111] directions are given below, 

7123(6/) = cxp[iJ'-6>] (i?31 J. + Rz2Jy + i?33^.) exp[-iJ'-0] 

(35a) 



/-I 



o\ 



jrwlll]) - 


1 

' 2 






-|V2^ -i 




Vo 


1/ 




/o 


V3 \ 


jr{[im - 


1 

' 2 





2 -^/i 

0-2 






lo 


-%/3 / 



(35b) 



(35c) 



Hence, the Zeeman interaction can always be written as, 

Hzeeman = -2^lB (« J3''' + 'Z^s'^)-^ (36) 



are com- 



sary for the description of the problem since 
plete. They are performed to obtain more convenient form of 
Hamiltonian. Examination of method of Eppenga and Schru- 
umans shows that they have simply performed the transforma- 
tion described by Eq.([3T]i and followed by summation over the 
Cartesian axis {1,2, 3}, but without the transformation of the 
bases to \(j)]^? j)- Since they are only dealing with symmetric 
Lowdin term without operator ordering, the result is linked to 
that of method of Luttinger by the unitary transformation of 



the bases describe by Eq.( 34 1 



VIII. COMPARISON OF DOUBLE GROUP 

FORMULATION WITH RELATIVISTICALLY 

CORRECTED DKK HAMILTONIAN 

Within the framework of one electron theory, the ap- 
proaches employed to solve the Schrodinger equation under 
the double or single group formulation should be equivalent 
(giving the same dispersion relations), provided the bases for 
perturbative expansion are complete. The primary difference 
between the two approaches can be identified as the inclusion 
or exclusion of H^^i in the unperturbed Hamiltonian. As a 
result, the bases in the perturbative expansion when solving 
Eq.lflll differ under the two approaches. 

In this section, we first examine relativistic corrections to 
the DKK model by method of invariants. The matrix repre- 
sentation i^loi and Hl^2 ^r^ then evaluated with respect to a 
complete set of single group product bases. The validity of 
DKK model with relativistic corrections are then examined in 
terms of its ability to explain experimental observations. 

The effects of spin orbit interaction can be incorporated into 
a limited band model using the method of invarianP and per- 
turbation theory '". Here we only consider relativistic correc- 
tions to the 6-band model. Under the method of invariant, 
possible perturbations due to the relativistic corrections are 
constructed from products of the matrix representation of the 
spin operator S, and components of wave vector fc. With re- 
spect to single group product bases, the matrix representation 
of spin operator S is simply |(T ® Inxn, where cr are Pauli 
matrices, and n is the dimension of the single group IR. Util- 
ising the character table of the Oh group, it can be shown that 
F+ ® F+ = {r+ e F+ ® F+} e [F+]. Hence, the symme- 
try allowed intra-band irreducible perturbations must have F^ 
symmetry, if it is time reversal odd, and F^, F;^, or F^ sym- 
metry if it is time reversal even. The permitted irreducible 
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TABLE V. Relativistic irreducible perturbations containing first order in Pauli matrices (a^), and zero or second order in wave vector compo- 
nents (fc^). Terms with FJ symmetry are time reversal odd, and terms with F^ , Ff, and F^ symmetry are time reversal even. 



IR 



Irreducible perturbation 



r+ 



ly^Bx + CTyBy + (JzBz 

J -^{2azBz ~ a^Bx — (TyBy) 



a.S, 


— (JyBy 


CTx 


oxk^ 


ay\ 


• ayk'^; 


f^z 


[azk' 


OyBz 


+ azBy 


^zBx 


+ cJxBz 


ax By 


+ (7yBx 



(Jyky , 



ayXkx: r^y\ ~r azxkz, rCxf 

^zxky, kz\ 'T (^xykx^ tZy\ 

^ ^x \kz^ kx / ~r <7y I ky , /Cz I 



perturbations, which are first order in a but zero or second 
order in fc, are listed in Table [V] 

The 6-band Hamiltonian with relativistic corrections con- 



sists of a non-relativistic term i/„r(fc), the k independent rela- 
tivistic correction Hr{0), and the k dependent relativistic cor- 
rection Hr{k), 



H{k) = Hnr{k) + HriO) + Hr{k) 
Hnrik) = Ul^ (12X2 ® HoKxik)) C/p+ 



r+ 



Hr{0) = - tf/^+ (a. ® 4 + c.p.) L/p+ = -f C/L (<T • I) C/r+ 

Hrik) = -^Ul^ (-aifc2 {a -I)- Sa^ [{a, ® h - a ■ I/i)kl + c.p] 

+ -a3[(o-:r (g) ly + cry(E) Ix){kx, ky} + c.p.] j C/p+ 
-HbUI+ {36i(a,B, + c.p.) (g) laxa + 8662 (((t,B,) g, [ll - 1^3) + c.p.) 

+663 {{(JxBy + <7yBx){Ix,Iy} + C.p.)} [/p+ 



(37a) 
(37b) 

(37c) 



(37d) 



where Hjjkk is given by Eq.(Dl 1, 11^+ defines the unitary pressed as 



transformation from single group product basis to adapted 
double group basis, c.p. denotes cyclic permutation of the pre- 
ceding term, a^ are the Pauli matrices, and I^ is the matrix 
representation of angular momentum (in units of K) with re- 
spect to the / = 2 real bases of F^ IR shown in appendix [aI 
They are given by. 



Ix^ 



'^0 0^ 

i 

\Q -i Oy 



/„ = 




71 = 7"'' + oi 


(39a) 


72 = 7r + 02 


(39b) 


73 == 7r + «3 


(39c) 


K = k"'' + bi + 762 - &3 


(39d) 


g = 9&2 - 363 


(39e) 


7r = ir " 2ai 


(39f) 


k'" = 4k"'' ^ 2bi - 3262 - I663 


(39g) 


72^^=72"'^-y 


(39h) 


73"^ = -73"'' + y 


(39i) 


^vs ^ ^nr ^ 25i - 462 - 263 


(39j) 



where the non-relativistic parameters are determined by the 
second order interaction parameters F, G, Hi and H2, as 



The relativistically corrected Luttinger invariants are then ex- shown in Eq.(D4i. Under this particular treatment of spin or- 
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bit interaction, it is clear that the existence of the q invariant 
is purely due to the relativistic correction. Moreover, the Lut- 
tinger invariants describing Lowdin interactions in different 
blocks, now differ from the non-relativistic values, and from 
each other. 



In general, most implementations of the single group for- 
mulation neglect the relativistic k dependent correction to the 
Luttinger invariants. Inclusion of Hr{k) in the Hamiltonian 
modifies the Luttinger parameters, and removes the link be- 
tween them in different blocks under the adapted double group 
bases, as originally envisaged by Kane^. It is therefore im- 
portant to establish the numerical impact of both Hr{k) and 
H,.{0). 



single group product bases of IRs a^, a" may be written as. 



(raA^k.p v^^) = E^^'f^"","" 






^l 



u^ 



I2 



M„ 



(40) 



where A/^u^i, is defined by Eq.(B2b 1 for the momentum op- 
erator. The matrix representation of i/fo2 ^i'^^ respect to the 
same bases may be written as. 



Under the single group formulation, in which the unperturbed 
Hamiltonian is Hq, the gradient of the potential may be ex- 
pressed as commutator of the momentum operator and Hq, 



Let us now examine the matrix representation of -fffoi, and 
Hi — Hf, + i?|o2 ^i'^h respect to the single group prod- 
uct bases. The matrix representation of Hf, with respect to 



VV{r) = j[p,H^]='-{pW,-H^Qp) 

Expressing the single group product bases in terms of their 
spinor and orbital part IV'q") — I '/'a") Ix*)^ ^ rn^Y be evalu- 
ated as. 



4rnoc2 



(C"|(^xVl/W) 



^i^ 



3 \ - 



V u „ = 



ix'l^W) ^ {€^\{pH^ ~ H^p)\<j>l.) 



ih 
AttIqC^ 

^- — 'a X M„..„. 

4mQC^ 



(41) 



r 



Hence, the matrix representation of Hi may be written as. 



<'"" = EfcA.(e 



Kt 



(42) 



Comparing Eq.(40i and Eq.(41 1, one can see that the re 



duced tensor matrix element of V and U differs by a factor 

S„.;-_E, 



^^-^ ^^, which is inversely proportional to the electron 



rest mass and is of the order of 10^^. This means that H^^2 is 
relativistically small compared to H^ , and may be neglected 
when the two terms occur together as part of Hi . 

The same can not be said about the k independent spin or- 
bit interaction H^^i. If all the zone centre states constitute a 

complete set of bases and i — 'J2af> n I'f'ai^) {'^^f> I' ^^^^ ^^ 
relevant matrix element with respect to the single group prod- 
uct bases is given by. 



(^^„|^.(VnOxp) 



e)- 



'imkc^ 



irrinc'^ 



{x'\^ W) ■ (e.|P^oi ^P-H^P^ X Pico 



imlc^ 



{x'\^W)-{ E (^UpK^) X {<t>Up\4>l^)E^, 



E('^'"ipic^)x(oipiC")^o'' 
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Then the matrix representation of H^^i with respect to the 
single group product bases, denoted as Wa^ ,a^ may be written 

as. 



W^. 



Am^c^ 



O- • ^ Pa-,af> X PaC.a- {Eaf " Ea^ 



Ah'^c^ 



a- ■ J2 M^-M' X -/^a^a"(^a^ - EiA}) 



where P is the matrix representation of momentum opera- 
tor The inverse rest mass dependence of the matrix represen- 
tation of Hl^i is absorbed into known, finite reduced tensor 
elements of M matrices. It is clear that this term is not rela- 
tivistically small from this particular view point. The angular 
dependence from cr ■ (Af^u.^^ x M^p.,^^^ is readily diago- 

nalised as a block, by a unitary transformation such as U^^ , 
given in appendix |C] It is important to recognise that while 
intra-band terms become diagonal as part of the Hamiltonian 
after the unitary transformation, the inter-band k indepen- 



dent spin orbit interaction terms remain on the off-diagonal 
blocks when evaluated with respect to the adapted double 
group bases. These inter-band terms are generally ignored 
in models derived under the single group formulation. Since 
Vy(r) X p transforms according to F^ IR of the Oh group, 
single group selection rules readily permits this k independent 
interaction as intra-band terms between states with F^ or F^ 
symmetries, or inter-band terms between states with F3 , F4 
or Fg symmetries with the same spatial parity. Hence, one 
can not have a finite reduced tensor elements for the intra- 
band terms, without having a finite reduced tensor elements 
for inter-band terms. 



The k dependent relativistic correction arises from con- 
struction of Lowdin terms using the first order k ■ n pertur- 
bation. Using the matrix representation of the Hi given in 
Eq.(42i in place of Af^u ^i. in Eq.( 17c|17d| to evaluate the 
Lowdin interaction with respect to the single group product 
bases, we obtain. 



x^y.z x,y^z 
ra",a'" __ ST^ ST^ ST^ 

Lowdin / j / j / j 

a/3 At V 



{Ea^ + Eav — 2Eaf>) 

4{Ea^ - E^,){E^. - E^,)' 



{k^^.^u:, 



± Ua^,af<U^p 



non-relativistic 



+ \^l^,aK^af>,a- + ^a-^ .aP^af ,a-) =•= (^a" .a/^ ^q/3 ,a- + ^a" .a^ f^Q/3 ,q- j 



relativistic correction; first order 



ya^.af<^af>,a^ ± KT" .a/^ V"^/j ,„„ 



relativistic correction: second order 



(44) 



r 



where U and V are those defined in Eq.(40 41 1, and the ± 
signs refers to symmetric and antisymmetric Lowdin interac- 
tion. It is apparent from Eq.(J4T]i that V is zero when a" = a" 
as demonstrated by Bir and Pikus'*. Hence, there is no intra- 
band contribution to the relativistic correction to Lowdin term. 
The first order relativistic corrections to the Lowdin interac- 



tions, shown in Eq.(44 1, have reduced tensor elements which 
are a factor of °im c^" smaller compared with the corre- 
sponding non-relativistic Lowdin interaction. This indicates 
that the k dependent relativistic correction to the Lowdin in- 
teraction can be neglected when treating the problem using 
single group product bases or adapted double group bases. 

While the method of invariant shows that the Luttinger in- 
variants differ between blocks in the 6 band model, the rela- 
tivistic corrections that lead to these differences are negligi- 
ble. Thus, the equality among these Luttinger invariants in 
different blocks remains valid. The q invariant is permitted to 
exist by symmetry, but is very close to zero. Therefore, the 
DKK model with relativistic corrections still can not explain 
phenomenon readily observed experimentally. For example. 



experiments indicate a non-zero q invariant, and an effective 
mass for the spin orbit band that differs from the value ob- 
tained from 71. In addition, the perturbation method could not 
account for the mixing introduced among states in the remote 
set by the k independent spin orbit interaction. In terms of k 
dependence, there are four independent second order interac- 
tion parameters in the valence band, but five Luttinger invari- 
ants obtained from method of invariant using double group 
selection rules. 

Within the single group formulation, it has been assumed 
that the k independent spin orbit interaction with respect to 
the single group product bases is finite for intra-band terms in 
order to explain the experimentally observed splitting (A) due 



to spin orbit interaction. The relation shown in Eq.(43 1 relates 
the reduced tensor elements of matrix representation of H^^^ 
to reduced tensor elements of momentum operators. Thus, a 
finite intra-band term would imply that the symmetry permit- 
ted inter-band term is also finite. It is then important to include 
such terms in any model. However, the inter-band terms are 
generally ignored in single group formulations. Furthermore, 
its effect could not be included in the treatment of spin orbit 
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interaction under limited set of zone centre states, as shown in 
the following paragraph. 

Inter-band k independent spin orbit interaction terms cause 
mixing between states, particularly those with F^ symmetry 
but derived from T^ and FJ single group states. Unless these 
terms are included explicitly in the single group formulation 
of multiband models, the implicit effect of mixing leads to de- 
parture of zone centre states from their single group heritage 
of one specific single group IR. In other words, it calls for 
the double group selection rules and formulation. If states 
with Fg^ symmetry are not included explicitly in the near 
set, then such an effect may manifest itself in the absence of 
equality between Luttinger invariants in different blocks ac- 
cording to double group classification but otherwise associ- 
ated via its single group heritage. In anycase, the universal 
constants at the front of the Eq.(|43]l still give a scaling fac- 



tor of 6.4 X lO^^eF^A . Thus, contributions from remote 
states beyond the normal anti-bonding zone centre states are 
required for the summation in Eq.([43]l to yield the level of 
spin splitting observed in experiments. Such remote states 
would have energies which invalidate the assumptions of the 
Foldy-Wouthuysen transform. This raises a question on the 
validity of the general assumption that the reduced tensor el- 



ement in Eq.(43 1 is finite for the single group product bases. 
Alternatively, one may question the validity of single group 
product bases as representations of the zone centre states. A 
mixed bases, wether caused by inter-band k independent spin 
orbit interactions or many electron effects in the construc- 
tion of zone centre states, would support a finite spin splitting 
and k dependent relativistic corrections. Such mixed bases 
would then requires the double group selection rules and dou- 
ble group formulation. 

The Kane model'"* ^^ and extended Kane model '^'^^' were 
developed to account for finite spin orbit splitting, and to in- 
corporate three or more zone centre energy levels in the multi- 
band model. It assumes the zone centre states in the near set 
have energies of those obtained experimentally, and are there- 
fore classified according to double group IRs. The model re- 
tains single group selection rules and reduced tensor elements 
of the single group first order interactions. Implicitly, no mix- 
ing is permitted upon the formation of zone centre states clas- 
sified according to double group IRs. The effect of the remote 
set of states is included by using experimentally measured 
Luttinger parameters, which are then modified by the removal 
of contribution from states newly included in the near set ac- 
cording to their double group classified energetic positions. 
The relations between Luttinger invariants shown in Eq.(|D4f| i 
are assumed to be true, but the relation shown in Eq.( |D4e i is 
not. While it is clear from the procedure of modifying Lut- 
tinger parameters for the 8- or 14-band models that the inter- 
band parameters would be different from the valence band due 
to different energies, this distinction is not made under either 
the Kane or extended Kane models. On the other hand, the 
distinction was made by Weiler et al.ES under the same single 
group selection rules, but with double group classification of 
remote as well as near set of states. The impact of this dis- 
tinction is explored in the following section by extending the 
Weiler model to include additional bands, and comparing the 



results with those of experiment. 



IX. COMPARISON WITH WEILER MODEL AND 30 
BAND MODEL 

In models with a limited number of bands, the magnitude of 
relativistic corrections to Luttinger invariants from the DKK 
model were found to be negligible. As a result, the relations 
between Luttinger invariants and second order interaction pa- 
rameters differ between the single group formulation with rel- 
ativistic corrections, and double group formulation. Specifi- 
cally the lack of a q invariant was problematic in the descrip- 
tion of valence band structure. The Weiler model'^ was intro- 
duced to account for some of these experimental facts. The 
Weiler model assumes that the zone centre states have ener- 
gies of those obtained experimentally, and are classified ac- 
cording to double group IRs. It however retains single group 
selection rules, and the reduced tensor elements of first order 
interaction in related blocks derived from that of the corre- 
sponding single group first order interaction. This treatment 
of the reduced tensor elements is identical to the approach 
employed in the extended Kane model for states in the near 
set. Thus, the first order interaction matrices are the same 
as those obtained by transforming the first order interaction 
under the single group into those under the adapted double 
group, as described in appendix IE] In this respect, the Weiler 
model is still based on single group formulation. The model 
is ideally suited to explain the finite q invariant in the valence 
band block through exactly the same mechanism proposed by 
Hensel and Suzuki"*^. However, the links between the double 
group classified zone centre states, and single group selection 
rules or reduced tensor elements of first order interactions, are 
ambiguous. While it recognises that Luttinger invariants in 
different blocks of Hamiltonian are different (as opposed to 



Eq.(D4fi), no attempt was made to define independent sec- 
ond order interaction parameters and relate them to the rele- 
vant Luttinger invariants. In the Weiler model, remote states 
derived from single group states with Fg^ symmetry were ex- 
cluded. An implicit assumption of the model is the absence of 
inter-band k independent spin orbit interaction terms, and its 
effect of mixing of zone centre states. 

In this section, the Weiler model is extended (referred to 
as extended Weiler model in the rest of the manuscript) to in- 
clude states with F^ (F^^ ) and Fg © F^ (F^ ) symmetry in the 
remote set of states. Relevant independent second order in- 
teraction parameters are defined taking into account the dou- 
ble group classification of zone centre states in terms of their 
energies, but retaining the single group first order interaction 
parameters which adhere to single group selection rules. Rela- 
tions are then derived between Luttinger invariants and second 
order interaction parameters so that material parameters can 
be extracted under the assumptions of Weiler model, and their 
values can be compared with those of the experimentally de- 
termined parameters. Issues relating to the assumptions made 
in the Weiler model and other single group formulations, are 
addressed in the next section. 

It is the case that many of the effects of the double group 
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FIG. 1. Schematics of ordering of energies of zone centre states, allowed first order interactions and their interaction parameters under (a) 
double group formulation and (b) 30 band single group formulation for Ge. 



formulation pertaining to the energetic positions of the zone 
centre states, are already incorporated into the Weiler model 
through the use of empirical fitting parameters. While the en- 
ergy of zone centre states with Tj symmetry are unaffected 
by spin orbit interaction in the single group formulation, the 
use of empirical parameters allows any such effect to be in- 
corporated in a practical implementation of the Weiler model. 
As shown in previous sections, the angular dependent form 
of the Hamiltonian for models up to, and including 8-bands 
are equivalent to those of the double group formulation. With 
parametric fitting, there then appears to be little room for dif- 
ferences to arise between the results of the double group for- 
mulation and the extended Weiler model. 

However, the main distinction between the two formula- 
tions is the allowance (or exclusion) of mixing between zone 
centre states derived from different single group symmetries. 
The effects of mixing in the single and double group formu- 
lations are exemplified in FigfTl where the energetic positions 
of zone centre states and allowed first order interactions are 
shown schematically. Under the single group formulation, 
there is a one-to-one correspondence between the single group 
zone centre bases and the adapted double group IR (or pairs of 
IRs where spin orbit splitting exist). Allowed interactions are 
those permitted by single group selection rales. There are two 
inter-band k independent spin orbit interaction terms identi- 
fied as Ap. and AfJ. The term Ai, in particular allows for 
mixing discussed at end of section III and their incorporation 
would ensure the equivalence between single group and dou- 
ble group formulations based on one electron theory. How- 
ever, this term is either set to zercPSlor not included at all-in 
the 30-band model, and is absent in any other implementation 
of single group formulation. Consequently, mixing caused by 
inter-band k independent spin orbit interaction is ignored. If 
this term is included, it would introduce approximately 4% of 
Fg (r^" ) state into the Fg" (F4 ) state, assuming a energy sepa- 



ration E^, - E^, ^ 7.27eV^'^ and A^^ = 0.29eF in Ge. The 
effect of the inter-band k independent spin orbit interaction 
amounts to lifting of restrictions in the adapted double group 
bases and necessitates the double group selection rules. 

In the double group formulation, the zone centre states are 
solutions of unperturbed Hamiltonian Ha in Eq.lfTli, which al- 
ready contains the k independent spin orbit interaction term. 
Therefore, the zone centre states are already mixed in terms 
of adapted double group bases. Consequently, first order in- 
teractions are those permitted under the double group selec- 
tion rules. Specifically, the interaction between states with 
Fg IRs is characterised by two interaction parameters and two 



linearly independent matrices as described in Eq.(B5 1. To ad- 
here to double group selection rules, the two linearly indepen- 
dent matrices which form the first order interaction must be 
combined when the Lowdin interaction is constructed using 



Eq.( 17d I. This introduces the three linearly independent sec- 
ond order interaction parameters, C + '„+ ' C+ V+ '^^'^ '^r+ 'r+ 



5pi- p-f 



± 



for Lowdin interactions between states with Fg symmetry. In 
contrast, any single group formulation only permits zone cen- 
tre states to be associated with one single group IR. In this 
case, the first order interaction between states of Fg IR can 
always be described by one interaction matrix if states derived 
from Ft^ , F;^ , F J symmetry are excluded. The contributions 
to the Lowdin term can be added at 2nd order level since each 
occurrence of these interactions is always associated with a 
different set of remote states. 

The angular dependent parts of the first order interaction 



associated with ^^^ 



and ^^_|_ double group parameters. 



are the same as the angular dependent parts of the interaction 
associated with P^ and P^ single group parameters, respec- 
tively. The distinction between single, and double group for- 
mulations is evident when constructing a Lowdin term. The 

_ for a given pair 



interaction described by ^^^ 



- and e+ r- 
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of Fg IR must be added at the first order level, whereas the 

corresponding terms PP and P^ must be added at the sec- 

r^ 3 
ond order level. Thus, there is no equivalent of the (l ' + 

parameter under the single group formulation, unless states 
with r^ © r^(r^) symmetry are included as remote states. 
The main difference between EWZ model and the extended 
Weiler model with 14-, or more bands is explicit in the num- 
ber of linearly independent matrices required to describe first 
order interactions between states with Fg symmetry ^. The 
differences for models with fewer bands are implicit in the re- 
lations between second order interaction parameters and Lut- 
tinger invariants. The main impacts are then to the relations 



between the Luttinger invariants in various inter-, and intra- 
band blocks in the Hamiltonian, the fitting of experimentally 
derived Luttinger parameters in these different blocks, and the 
dispersion relations away from the zone centre. 

Under the single group selection rules, the first order in- 
teractions between states are given by Eq.( |B3[ ). The corre- 
sponding first order interaction between adapted double group 
bases may then be obtained using the procedure described in 
appendix IE] While it is not possible to determine the sign of 
the first order interaction parameters (^) from symmetry argu- 
ment, the sign of ^s between states with double IRs derived 
from the same parent single group IR, are related. The fol- 
lowing relations can then be extracted. 



Pp -Cr+(p),r-(s* 
Pp =^r+(p),r-(p* 



Pf = ^Tf 



jd* 



Pp =^r+ 



r;r(rt^r-(/*) 



(p),r3-(<i*) 



= "^6er+(r+),r-(r-) = ^/3e^+(^+),^-(^-) 
= -v3Cr+(r+),r-(r-) = ^^^^r-(r-),r+(r+ 
= ^^r+(r+),r~(r-) = ^^r-(r-),r+(r+) 
= V6^r+(r+),r-(r-) = ^^r-(r-),r+(r+) 
= -v6Cr+(r+),r-7(r-) = v3er+(r+),r~(r-) 
= ^^r+(r+),r-(r-) = ^^r+(r+)T3-(r3-) 



(45a) 



(45b) 
(45c) 
(45d) 



The interaction between states with F^ symmetry with respect to the adapted double group bases may be written in block form 



as. 
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2x2 



'K^t,r-)Ur 



\n 



r+,r, 
K 



K.I 



r7'+,r: 



-) K. 



r+,r- 
^r+.r; 



(46) 



with relations between first order parameters given as. 



-> ^^r+(p)T-(/*) 



'^r+(r+),r-(r-) 



^^r+(r+),r-(r3-) 



r8+(r5+).rf(r5") 



-3\/2^, 



r?(r5^)A"(r5") ~ 2^r?(r+)T7(r-) 



?P+^T^+^ 



(47a) 
(47b) 



r 



Since the first order interaction parameters ^ are all rear^, 
these relations permit some of the sign differences between re- 
lated parameters to be fixed. These relations shall prove useful 
in determining the sign of second order interaction parameters 



in the inter-band block shown in Eq.(48 i, within the extended 



Weiler model. 

Having identified all the possible first order interactions 
under single group selection rules, second order interaction 
parameters due to Lowdin interaction can be defined in the 
same way C's are under double group formulation of EWZ, or 
F, G, Hi,H2 are under the single group formulation of DKK. 
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(48a) 



(48b) 



(48c) 



(48d) 



(48e) 



(48f) 



(48g) 



(48h) 



(48i) 



(48j) 
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where the single group interaction parameters are defined in 
Eq.(45|46i. The appropriate negative sign in front of terms 



p^ p^ 1 

such as Z I + and Z + ' + are derived from the signs given 

^ S '^ 7 ^ 8 '^ 7 



inEq.(|45i. 

With these definitions, symmetric and anti-symmetric 



Lowdin terms are constructed using Eq.(17i. The resuhing 
Luttinger parameters in the 6 band model may be expressed 
in terms of these 10 independent second order interaction pa- 
rameters as follows. 
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(49a) 
(49b) 
(49c) 
(49d) 

(49e) 

(49f) 
(49g) 
(49h) 

(49i) 

(49j) 



In addition, the conduction band Luttinger invariants may be 



written as 



7r 



zl 



T- .V- 



4z: 



V- .V- 



2Zl 



4z: 



50 



(49k) 
(491) 



In Eq.(48 i, contributions from remote states with Fg sym- 



metry are included in curly braces. The first order interac- 



tion /C 



r?(rj)^r-(r-) 



-,p-s is expressed as an linear combination of 



/Cj,+ and /C?,+ , as shown in Eq.(|46|l. In this respect, we 

^ 8 '^ S ^ 8 '^ 8 



may label states of Fg symmetry derived from Fg as 'mixed' 
with respect to the bases of Fg^ derived from r4 and Fg . 

p^ -3 

Such contributions add the additional parameter of Z \ _,_ 

'-&''-& 
to the valence band block. It is clear that these parameters 

have the same labelling scheme as the C, parameters defined in 

EWZ, and there is a one-to-one correspondence provided that 

remote states with Fg^ symmetry are included. Thus, the pres- 



ence of Z, 



r7,3 



in the extended Weiler model is dependent 

Simi- 



-p+ p+ 

on the inclusion of remote states with Fg" symmetry. 
lar arguments can be made about single group states with F;^ 
and F^ symmetry which give rise to double group bases of 
F^. In the implementation of Weiler model or 30 band model 
in the literature, remote states with Fg" symmetry are always 
excluded. Hence, there are only four independent second or- 
der interaction parameters in the valence band related to the 
five Luttinger invariants. As one expects the invariants to be 



r 



linearly independent of each other, the presence of only four 
parameters in Eq.(49i remove such linear independence. In 



addition, the Weiler model has difficulty in explaining the ef- 
fective mass in the spin orbit band of Si and conduction band 
of Ge as shown below. 



The relations shown in Eq.(49l derived from the extended 



Weiler model may be compared with Eq.( |Dl|D"4l l derived 
from the DKK model with relativistic corrections, and those in 



Eq.(23 I derived from the double group formulation. The Lut- 
tinger invariants in each of the blocks are now different from 
each other The magnitude of these difference is determined 
by the different energy scaling in the corresponding Z param- 
eters. The q parameter is now non-zero, provided spin orbit 
splitting exists in the conduction band state derived from the 
antibonding p* state. It is therefore clear that extended Weiler 
model overcomes some of the key difficulties encountered in 
the DKK model with relativistic corrections. 

The impacts of classifying the energetic positions of the 
zone centre states according to double group under the Weiler 
model, is best seen in the contraction of the full zone 30-band 
model^**^^ into models with fewer bands. The relations be- 
tween the 30-band single group interaction parameters and 
Luttinger invariants, have been obtained by Richards et al.^^ 
and are given here for crystals with diamond lattice. These 
relations have been obtained by transforming the Hamiltonian 
referring to product bases^'-* into adapted double group bases, 
using the diagonalization matrices shown in appendix lO The 
Lowdin terms are then constructed using the resulting interac- 



tion matrices and parameters, and Eq.( 17 1 giving the follow- 
ing relations. 
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(50a) 



(50b) 



(50c) 



(50d) 



(50e) 



where the definition of the matrix elements (reduced tensor For the spin split off band under 2-band or 6-band model. 



d* 



element) are shown in FigjTp, and go ~ 2. The effect of A 
is not included, in accordance with the single group formu- 
lation in the literature. The original atomic orbital has been 
shown as superscript in the energies of the relevant zone cen- 
tre states. The order of basis for valence band used in the 
original work of Richard^^, is the same as those in Eq.( 15 1, 



and hence there is a difference in sign of the 72 invariant. The 
contribution from each remote state symmetry is described by 
Winkler^. The number of parameters, including the zone cen- 
tre state energies, are large, but they can be reduced to the set 



of independent Z parameters defined in Eq.(48 1 



Referring to the 30-band model and interactions designated 
in Figure [T[ the Z parameters can be related to single group 
first order interaction parameters as follows. For the valence 
band under 4-band model. 




(5 If) 



(51g) 



For the conduction band with T^ symmetry under 2-band 
model. 
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For the inter band block under 6-band model, 
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There are no states derived from Fg symmetry in the 30 band 

model and therefore the parameter Z + ' , is absent. 

The implementation of the 30-band modeFS122l was in- 
(51e) tended to provide both an accurate description of dispersion 
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over the whole of the Brillouin zone, and the extraction of 
energy gaps at the other high symmetry points. It is how- 
ever the case that the Luttinger parameters, which are derived 



from the first order parameters in Ref58 using Eq.(50i, de 



viate substantially from the experimentally measured param- 
eters. This means that the dispersion relations calculated near 
the zone centre are not accurate. Utilising the Weiler model. 
It is possible to extract a set of single group first order interac- 
tion parameters from the experimentally determined Luttinger 
parameters of the valence band. This then permits the extrac- 
tion of a full set of Luttinger parameters for the inter-band, 
and spin split-off band blocks, and places a lower limit to the 
Luttinger parameter in the conduction band (Fy ) block. These 
derived parameters can be compared with experimentally de- 
termined values, where available, to ascertain the veracity of 
assumptions in the Weiler model. 

The single group first order interaction parameters can be 
extracted from experimentally determined Luttinger parame- 
ters 71,72 and 73 with the aid of Eq.( 51|49 l and some ad- 
ditional assumptions. Ignoring the remote anti-bonding /* 
state with r^(r2^) symmetry''', the first order interaction pa- 
rameters P;* (er+ r- ) > Pr (Cr+ p- ) and Pi** {^j.+ p- ) can be 



obtained from Eq.(5Tj49]) given the zone centre state energies. 
Then a full set of Luttinger parameters in the inter-band block, 
spin split off band block can be obtained. The lower limit to 
Luttinger parameter in the conduction band (Fy ) block can 
also be evaluated by ignoring the contribution of d states with 
F^ © F^ symmetry in the conduction band. These parameters 
are listed in TablelVIIfor Si and Ge. 



TABLE VI. Luttinger parameters of different blocks derived from 
71 , 72 , 73 of the valence band under the Weiler model 





71 


72 


73 


72 


73' 


7r 


7?(rf) 


Si 
Ge 


5.26 
13.38 


-0.38 
-4.24 


1.56 
5.69 


-0.39 
-3.73 


-1.55 
-4.98 


5.18 
10.47 


> -5.76 
> -21.91 



Under the 2 band model, the effective mass of the spin or- 
bit band is smaller compared with experimentally measured 
value in Si, and the lower limit of conduction band effective 
mass is larger compared with experimentally measured value 
in Ge. In addition, the dispersion along the (110) direction 
in Si, under the 6-band model is un-physical (see Figure l4]l. 
These discrepancies are generally attributed to remote states 
that have not been included in the model. For example, the in- 
troduction of the Hermann- WeisbucH^ parameters accounts 
for a similar over estimation of the conduction band effective 
mass in compound semiconductors. It is however clear, that 
all of the valence band states are well identified by experimen- 
tal methods, and any excluded remote states reside above the 
lowest conduction band. Such states would include those in 
the conduction band with Fy (F^) symmetry from / atomic 
orbital, and states with F^ ® Fy" symmetry from d atomic 
orbital. The contribution from such un-known, or excluded 
states, to the Lowdin interaction would only increase the ef- 
fective mass of the lowest conduction band, and reduce the 



effective mass of the spin split off band. These discrepancies 
are thus intrinsic to the Weiler model. 

Under the extended Weiler model, it is possible to provide a 
set of second order interaction parameters to predict the appro- 
priate effective mass of conduction band in Ge and spin-split- 
off band in Si. The extended Weiler model includes remote 
states with Fg" Fy (F^r) symmetry, and the second order 



interaction parameters include Z 






However, the mag- 



r^ 3 r^ 2 

nitude of Z l ' + would exceed that of Z I ' + (see section 



XI I, which would not be consistent with known energies and 
symmetry of the zone centre states. 



X. LOWDIN INTERACTION IN MULTIBAND MODELS 

In both the double group formulation of EWZ and single 
group formulation in the general literature, the Lowdin inter- 
action terms can be constructed from first order interaction 
matrices and have exactly the same form. This is the case 
even though the basis functions for the perturbative expansion 
differ This section examines the way in which the Lowdin 
interaction is derived, and shows that it is only permissible to 
use the double group formulation to derive the solution of one 
electron Schrodinger equation shown in Eq.(fT]l when a limited 
set of bases is stipulated. 

It is well known that the k independent spin orbit inter- 
action term, i?|oi, has non-zero matrix elements as an intra- 
band term. Within the single group formulation, such ma- 
trix elements are also permitted as inter-band terms between 
states with F^, FJ, or F^ symmetry with the same spatial 
parity. With respect to the single group product basis, the ma- 



trix representation of this interaction was obtained in Eq.(43 1 
and shown not to be relativistically small. The Hamiltonian in 
Eq.(fT]l can be divided into three terms, 

with respect to a complete set of bases. To obtain the Hamilto- 
nian with limited multiband model (e.g. 8-band), the complete 
set of zone centre states are divided into near and remote sets, 
which have relatively weak interactions between them. The 
interaction between states in the near set are included explic- 
itly, and the effects of remote set on the near set are included 
as perturbation to the desired order. The Lowdin method calls 
for a unitary transformation which leaves the inter-band block 
between the near and remote sets to be zero within the desired 
order in the relevant perturbation"* *■. 

Here, we wish to include all the terms quadratic in k com- 
ponents. To achieve this, the resulting bases must diagonalise 
Ha = Hq + Hl^i{cr), rather than just H^. In absence of the 
inter-band k independent spin orbit interaction terms, there is 
no mixing between states derived from different single group 
symmetry. The unitary transformation described in appendix 
[C] which is frequently used in the literature, would then diag- 
onalise Hq and H^gi{a) with respect to the adapted double 
group bases. This is clearly not true if inter-band k indepen- 
dent spin orbit interaction is present. Such interactions can not 
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be made diagonal in the Hamiltonian by the unitary transfor- 
mation within each IR. DiagonaHzation of Hq leads to mixing 
between single group zone centre states which necessitate the 
double group selection rules and formulation. In this case, the 
bases for the Lowdin term diagonalise Hq = Hq + H^^i{(t) 
and differ from the adapted double group bases. They are 
the bases relevant to the double group formulation. Hence, 
the Lowdin interaction obtained from perturbation theory in 
EWZ and any method of invariant*^ obtained using any form 
of double group bases^-' actually refers to the bases of double 
group formulation in the presence of inter-band k independent 
spin orbit interaction. 

The relativistic corrections to the DKK Hamiltonian, dis- 



cussed in the section VIII refers to the adapted double group 
bases. The distinction between the relativistically corrected 
single group formulation and the double group formulation 
lies in the zone centre bases of both the near set and remote 
set of states. The effect of mixing due to k independent inter- 
band spin orbit interaction can not be accounted using the per- 
turbation method employed in section [VlII] unless such inter- 
action is between states of the near set. Under the 30-band 
model, the relevant inter-band k independent spin orbit inter- 
action can be identified as the Ai term (which correspond 
to A3C and is set to zero in Ref^F"). Terms such as A^, also 
cause inter-band mixing but do not change the symmetry of 
the bases. Hence, it has no direct effect on the form of first 
order interaction matrices. The limitation of the relativistic 
corrected DKK model, based on adapted double group bases, 
is therefore attributable to the failure in dealing with the effect 
of inter-band k independent spin orbit interactions. 

In the extended Weiler model, the ambiguity between the 
double group classification of energetic position of zone cen- 
tre states and the use of single group selection rule and single 
group derived interactions, means that the zone centre states 
are not well specified. In essence, it lacks a clearly defined 
un-perturbed HamiltoniarPS The use of experimentally ob- 
served zone centre states under a double group classification 
means the model closely resembles the double group formu- 
lation, as evident in the identical form of Eq.([49| and Eq.([23]l. 
However, the 'mixing' of states with the Fg symmetry is in- 
corporated in the extended Weiler model only when remote 
states with T^ symmetry are included. If these remote states 
with Fg^ symmetry are considered, parametric fitting of ex- 
perimentally measured Luttinger parameters should return the 
same results as the double group formulation. Nevertheless, 
the Z parameters may not make physical sense from pertur- 
bation theory perspective^^. However, remote states with F^" 
symmetry are normally excluded under the Weiler model or 
30 band model. In this case, there are only four indepen- 
dent second order interaction parameters for the valence band 

p— p— p— 1 p- i-f 

r+ r+^ 



ZJ^. „i , Z'Ji '^4- and Zli ''j^) but five Luttinger in- 



(■^P+ p+ : ""p+ p+ 1 '"p+ r + 
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variants. The mismatch between the four independent sec- 
ond order interaction parameters and the Luttinger invariants, 
leads to linear dependence among the latter in contradiction to 
symmetry argument. 

The existence of a finite inter-band k independent spin orbit 
interaction, means that they must be treated under the unper- 



turbed system if quasi-degenerate perturbation theory is to be 
used to incorporate effects of remote states in a limited multi- 
band model. Hence the double group formulation, in which 
the mixing is automatically permitted under double group se- 
lection rules and material parameters, should be used. The 
additional mixing mechanism arising from formation of hy- 
bridised orbitals with consideration of spin means that the use 
of double group bases is essential in order to obtain the correct 
description of electronic dispersion under the k • p method. 
The Hamiltonian obtained from method of invariant does not 
specify the bases used for k expansion other than its transfor- 
mation properties in terms of IR of the symmetry group and 
the order of the basis function. However, the k dependence 
is consistent with that of the first order k ■ tv perturbation and 
Lowdin interactions under the quasi-degenerate perturbation 
theory. As a result, the basis must be those of the double group 
formulation instead of the adapted double group bases. 

It is clear from the discussion that the decoupling of states 
in near and remote set in a limited bases set model should 
be performed in such a way that all k independent perturba- 
tion is diagonalised by the unitary transformation. Apart from 
spin orbit interaction, strain is another perturbation that is fe 
independent. In principle, the method of dealing with strain 
perturbation should be different from those prescribed in the 
literature^. Upon the application of homogeneous stress, the 
symmetry of the crystal is lowered. The appropriate selec- 
tion rules to treat the k • n perturbation should be those of 
the prevailing sub-group H derived from the point group G of 
the unstrained crystal. The IR of the point group G then gen- 
erally becomes reducible under the compatibility relation be- 
tween H and G. Issues may arise from 'mixing' caused by the 
strain perturbation, as this modifies the first order interaction 
parameter or Luttinger invariants derived from the unstrained 
crystal. One can see some analogy between mixing due to k 
independent inter-band spin orbit interaction with respect to 
the the adapted double group bases and mixing due to strain 
perturbation. The treatment of stress as described by Bir and 
Pikus'^ may be justifiable if the terms in the Hamiltonian due 
to strain and causing mixing is smaller compared to other hav- 
ing similar mixing effects. However, the mixing would not be 
negligible in the infinite stress limit employed by Hensel and 
Suzuki"^ ^* when extracting the Luttinger parameters from ex- 
perimental data. The assumption holds for stress applied the 
the (001) direction, since all the terms first order in strain are 
diagonal, but fails in the (111) and (110) directions. This may 
cast some doubt on the particular method used. 



XI. MATERIAL PARAMETERS FOR SI AND GE IN 
DOUBLE GROUP FORMULATION 

The single group formulation provides an efficient way 
of calculating dispersion relations, but has some deficiencies 
which were identified and discussed in the last few sections. 
These deficiencies are generally ignored in the literature but 
can be overcome by using the double group formulation. To 
implement double group formulation with a limited set of near 
states, both linear and quadratic k terms must be considered 
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FIG. 2. Parametric dependence of dimensionless double group material parameters of the valence band, Luttinger parameters ^,5,72"', 

-4.24, and 73 — 5.69 for Ge. The symbols associated with the parameters are 



and 73" on (^ ^ ^ compliant to 71 — 13.38, 72 

dl „+ (•), CpT p+ (■). Cp+ ''+ (A), dl '"+ (T), C^l '1+ (m), 7r (o), 7r(n), contribution to k{a) from anti- symmetric Lowdin interaction, 

Ig.ig ^8 '8 8 '8 8 '8 8 "8 

p— 3 

and contribution to g(V) from anti-symmetric Lowdin interaction. Panel (a) and (b) are for (^ ^ ' ^ > and panels (c) and (d) are for 

•"8 '^8 
p~ 3 

(" J ' ^ < 0. See text for methods used to obtain these graphs. 



in the construction of the k-dependent Hamiltonian. These 
interaction terms should be defined by first order interaction 
parameters appropriate for the zone centre states included in 
the near set, and Luttinger invariants that describe interactions 
mediated by remote states. Given the one-to-one relations be- 
tween the Luttinger invariants and corresponding second or- 



and derived under a single group formulation. As the single 



der interaction parameters, as defined in Eq.(23 1, the Lowdin 



contribution to the k-dependent Hamiltonian may also be de- 
scribed in terms of the second order interaction parameters (. 

The extraction of Luttinger parameters from experimental 
data generally relies on magneto-optical measurements with 
or without an externally applied stress. The works of Pidgeon 
and Brown''^, Suzuki and Hensel'"^'*, and Trebin et al.'^^ 
provide the methodology and description of such measure- 
ments. Analysis of data relies on the use of k ■ p models with 
fitting parameters that are specified under perturbation theory 



group formulation utilises relations given in Eq.(D4f 1, the ex- 
tracted Luttinger parameters are model dependent. Moreover, 
the linear dependence between the Luttinger invariants under 
the single group formulation is problematic, as there are in- 
sufficient second order interaction parameters to describe all 
of the interactions. In this section, we examine the experi- 
mentally derived Luttinger parameters for both Ge and Si, and 
obtain a set of first order (^) and second order (() parameters 
under the double group fromulation. Data is available for the 
valence band (T'^), and spin split-off band (Fy^) for both Si 
and Ge, and the lowest conduction band (Fy ) for Ge. Second 
order interaction parameters (Q are extracted from Eq.(23i, 
and dispersion relations are calculated under the 6-, 8-, 14- 
band models. The dispersion relation should be physical, and 



the derived parameters should satisfy Eq.(23 1, with the appro 
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priate second order parameters remaining positive as required 
by their definition. 

The first concern is the second order interaction parame- 
ters, derivable from Luttinger parameters, that describe the 
main valence band block. The Luttinger parameters'^ for 
Ge are; 71 = 13.35, 72 = -4.25, 73 = 5.69, k = 3.59, 
and q = O.Ofr^. In principle, the second order parameters 



can be obtained by solving Eq.(23 1, using all the experimen- 
tally determined Luttinger parameters. However, the ( pa- 
rameters obtained from these values do not have the appro- 
priate signs as required by their definition. Assuming the 
Luttinger parameters associated with the symmetric Lowdin 
term are more reliable, it is possible to obtain the paramet- 
ric dependence of the ( parameters with the assumption that 
there are only 16 zone centre states corresponding to T^, 
r^, r^, Tg, r^ and Fg" IRs. Thus, the C parameters de- 
scribing the interaction with the F^ remote states have the 



relation C + ' + 



= ± W C+ ' + Cr+ 'p+ ■ The ± sign depends 

8 '^ 8 ^ 8 '^ 8 

on whether the first order interaction parameters ^^+ and 

^^+ , have the same or opposite sign. Choosing (l . as 
the parameter to vary, two sets of parametric relations can be 

T~ 3 

obtained depending on the choice of sign of C ? ' + . A choice 



on the sign of C 



rj,3 

r+ r+ 



and values for the second order param- 



eters is deferred, until discussions are made on the parameters 
concerning the spin split-off, and inter-band blocks. These 
relations are shown in Figurel2|a) and (c), and calculated val- 
ues of K and q are shown in Figure l2tb) and (d). It is clear 
that the values of k and q obtained in the literature do not 
provide a consistent solution to Eq.(23i. The importance of 



double group formulation can be seen in the larger magnitude 

r^ 3 r~ 2 

of CI ' + compared with CI ' . . Under the extended Weiler 

1 S '^ 8 ^ 8 '^ 8 

r^ 3 
model, it would be difficult to envisage such a large Z I ' . 

i 8 ^^ 8 



compared to Z 



r-,2 

r+ r+ 



given the remoteness of states with Fg 



symmetry. 

The Luttinger invariants in the inter-band block can not be 
measured directly in experiment. While it is recognised that 
such parameters should be different from those in the main 
valence band block'^, all previous models employed in the 
extraction of Luttinger parameters makes the assumption of 
Eq.(D4f I. To complete the 6-, and 8-band models, a system- 



atic way of obtaining these parameters under the double group 
formulation is required. A 14-band model is constructed for 
calculating the bulk dispersion using double group selection 



rules and first order interaction matrices given in Eq.(B4i, 
with the constraint that the Luttinger parameters, obtained by 
contracting a such model, should yield the experimentally de- 
termined values of 71,72, and 73. In other words, the sec- 
ond order interaction parameters describing the valence band 
block should follow those prescribed by the parametric rela- 
tions described in the previous paragraph and shown in Fig- 
urd2](a) or (c). The Luttinger parameters in the lower band 
models can then be evaluated by appropriate calculation of 
the Lowdin terms from the 14-band model, or more simply 



by using Eq.(23 1 and removing relevant terms corresponding 
to states included in the near set. This model diff ers from 
the corresponding single group formulation'^ '^ ^^ ^° model in 
terms of the description of the first order interaction between 
states with Fg symmetry, the order of the basis in F^ IR, and 
the absence of k independent inter-band spin orbit interaction. 
The Hamiltonian, in absence of other external perturbations, 
is then given by. 



H{k) = Ho + Hk.^k) + HLodwinik) 



Hk-K 



^ Lddwin\^ ) — 
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(52) 



(53) 



/ 



(54) 



where Hq is diagonal and contains the energies of the zone 
centre states, and the first order interaction matrices are 
given in Eq.(B4i. The matrix is ordered in states with 
Fo,F,- ^- ^^ 



column numbers. The lowest valence band zone centre states 
with Tq symmetry is considered as a remote state. The 
Lowdin term only exists between the conduction band states 



6 '^ 7 i-"^ 8 



and F7 symmetry with increasing row and with Fg and Fg symmetry due to double group selection 
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^r+ r+ 



TABLE VII. Double group parameters for valence band of silicon and germanium 



c 



r 2 

^r+ r+ 



c 



^r+ r+ 



c 



r r 1 

^« '^7 ^ « '^ 7 






7i 7i (1 7 ) 72 



73 



Ge 

Si 



0.48 
0.375 



5.715 
1.500 



3.980 
3.00 



0.253 
0.157 



-1.002 
-0.687 



2.456 
0.0288 



2.268 
1.105 



-3.786 
-0.208 



-3.007 
-1.821 



0.757 
0.417 



10.53 
3.448 



-26.32 
-7.028 



-2.693 
-0.039 



-3.732 
-0.845 



rules. However, the inclusion of Lowdin terms, that arise 



Under the Weiler model (which exclude the remote states 



from interactions with the remote state with Fg symmetry, with Fg symmetry), the lower limit of Fy conduction band 



have no impact on the dispersion of lowest conduction band, 
valence band, and spin spUt off band. As we are primarily in- 
terested in the dispersion of these bands, no Lowdin terms are 

considered in the 14-band model. Several relations can now 

p— p+ 

be derived from these assertions, CI + = —C - -, and 

^p+ p+ ^p p ' 

^ 8 i'- 8 J- 7 '^ 7 

Utilising the relations shown in Eq.(23 i 



Sp+ p+ - 
we have. 



-c 



rr,r; 



^r7- 



Ti 



'^r+,r+ ^ 4 l^i" 



l + 4fe 



^r+ r+ 



(55a) 
(55b) 



There are now only three adjustable parameters in the model. 



Cp+ p+: sr+.r^' V 



rr,r: 



given that the energies of the zone 



centre states are found to be E-p- 



E^ 



= 0.887eV^, 



Er, 



Ep+ = 3.006eV^, Ej.- 

1- ft J- Q 



E^+ = 3.206eF, and 



E-p+~E-p+ — — 0.297ey. Themagnitudeof the first order in- 

^7^8 

teraction parameters, such as ^p+ p- may be calculated from 

8 ■ 7 

the corresponding second order interaction parameter such as 



c 



riTl 



If one choose the condition that C 



rj.n 



> 0, it can 



be easily verified that the dispersion in the (110) direction is 

un-physical (valence band maximum not at zone centre) for 

p— 
all the possible range of C+ ^+- Thus, it is the case that 



5p+ p-t 



^rr,3 



'r+ r+ 

'- 8 ''- 8 



+ have opposite 

1 8 i'- 8 



Cl ''^+ < 0, and consequently ^p+ p- , C^ 

signs, and therefore (^ +' +, and Cp+ r+ also have opposite 



signs. Using the relations shown in Eq.(45 i as a guide, one 



p^ r^ 1 

can further assume that C T + < and C^ ' , 



<0. 



There is still a range of allowed values obtained from 

p^ 
the adjustable parameters of C Jf. ^, ^p+ p-, ^p+ p- that 

Fg .Fg 7 ' 7 7 ' 8 

provide physical dispersion relations in the valence band. 
The magneto-reflectance measurements of Aggarwal^' pro- 
vide values of Luttinger parameters 7f(F^) = — 26.32 (toc = 
0.038mo), and 7f" = 10.52 (m^^ = 0.095mo). These ad- 
ditional constraints reduce the number of adjustable parame- 

p^ 
ters to just one-^ "f^ +. The resulting dispersion in the X, L, 

n 

and K directions are shown in Figure pi with the correspond- 
ing second order parameters given in Table VII Given the 



way the fitting parameters are chosen, the dispersion at zone 
centre necessarily reflects the experimentally determined Lut- 
tinger parameters for the valence, conduction and spin-split 
off bands. 



effective mass can be extracted solely from the Luttinger pa 



rameters of the Fg valence band using relations in Eq.(49i 



This still exceed the experimentally measured value of Ag- 
garwal, and requires the inclusion of the Herman-Weisbuch 
parameter Under double group formulation the correct F^ 
conduction band effective mass can be extracted from the Lut- 
tinger parameters of the T^ valence band using only rela- 



tions in Eq.(23 i. The reasons for this can be demonstrated by 
comparing the relations between the 71 and second order in- 



teraction parameters given in Eq.(23 1, Eq.(49 1, and Eq.(D4i 
For the Weiler model, the Z 



P+ p+ 

J- 8 '^ 8 



0. Thus C 



is the 



term that differentiate the double group formulation of EWZ, 
Weiler model, and relativistically corrected DKK model. A 
smaller conduction band effective mass requires a larger mag- 



nitude of C - - 

^F, ,F, 



-c 






-zt 






or (7 (in 



F 3 

Eq.(|D4|). Because of the negative contribution of ( ^ ' ^ 

^8 '^8 



to 



71 in Eq.(23l, a larger magnitude of C 



rlT2 



is possible un- 



der EWZ compared to Z 



pi- p^ 



or a in Eq.(D4i for a given 



.F+,3 



value of 71 . Hence, the mixing term (I ' + is the key to the 

'^ 8 •'^ 8 

success of the double group formulation of EWZ in correctly 
evaluating the interaction between the zone centre valence and 
conduction band states. Similar effects means that inter-band 
Luttinger parameter 72** and 7!** under EWZ are smaller than 
their counter parts in the valence band. This allows the dis- 
persion in the (110) direction to remain physical. The impor- 
tance of zone centre states in the conduction band with Fg^ 
symmetry is further demonstrated by the differences observed 
between the 6-, 8-, and 14-band models. Thus, the disper- 
sion in the valence band and conduction band away from zone 
centre obtained in 6-, or 8-band models differ from that of the 
14-band model. 

For Si, there is no experimentally determined value for the 
F^ conduction band effective mass. However, the effective 
mass of the spin split-off band is knowrP^lto be rus — 0.29mo. 
Parameters may still be extracted for Si using a similar ap- 
proach to the one employed to extract parameters in Ge. A 
suitable set of material parameters for Si are shown in Table 
|VII| The dispersion relation for the valence and spin split off 
bands are calculated using double group parameters under 6-, 
and 14-band models, and compared with results from param- 
eters obtained under the 6-band Weiler model (Table |VI| l in 
Figure [4] The parameters shown in Table VI are used in the 
6-band Weiler model. The dispersion relation in the (110) di- 
rection shows the inadequacy of the single group formulation. 
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FIG. 3. Dispersion in Ge calculated using double group formulation, for 14-band (solid line), 8-band (dash line), and 6-band (dot-dash line) 
models. 



both in terms of the position of valence band maximum, and 
value of zone centre effective mass for the spin split-off band. 



The use of the double group formulation of the k-p method 
can provide a good agreement to a range of experimentally 
measured effective masses and thus provide an adequate de- 
scription of dispersion. However, one must be aware that ex- 
perimentally extracted Luttinger parameters are obtained us- 
ing models based on the single group formulation. Since all 
the valence band parameters for Ge were obtained using single 
group formulation'''*, there is no reason to give preference to 
the accuracy of Luttinger parameters associated with symmet- 
ric Lowdin term over those associated with anti-symmetric 
Lowdin terms. The discrepancy between the predicted Lut- 
tinger invariants (k and q) associated with anti-symmetric 



Lowdin term using Eq.(23 i from their experimentally deter- 
mined values is then due to the use of single group formulation 
in the analysis of data and uncertainty in the contribution from 
spin degree of freedom for mixed zone centre states. Never- 
theless, this is a significant improvement on the Weiler model, 
which has limited parameter set, and has difficulty in explain- 
ing effective masses measured in different bands simultane- 
ously. 



XII. DISCUSSION 



The treatment of spin orbit interaction is central to the dif- 
ferences between the double and single group formulations of 
the k • p method within the framework of one electron theory. 
The two formulations differ in the selection rules applied to 
the construction of matrix representation of fc • tt perturbation. 
They are equivalent when the zone centre bases for perturba- 
tive expansion in k are complete (no Lowdin interaction), as 
the only difference is in the choice of bases that describe the 
system. In this case, the spin orbit interaction terms, H^^^ and 
i7|o2' appear explicitly in the Hamiltonian. In constructing 
the Hamiltonian with a limited bases set, unitary transforma- 
tions are applied to all zone centre bases to eliminate the k ■ n 
perturbation between states in the near and remote set to the 
desired order, and diagonalise the un-perturbed Hamiltonian 
if ■ The double group formulation of EWZ treats ifo as its 
un-perturbed Hamiltonian, with the effects of H^^i incorpo- 
rated within the zone centre states and reflected by their ener- 
gies. Since these eigenstates form the bases of representation 
of the double group, double group selection rules are required 
in the description of fc • tt interaction between zone centre 
states. In contrast, the single group formulation retains the 
single group selection rules when evaluating the k ■ n pertur- 
bation. As such, the adapted double group bases are derived 
from unitary transformations of bases within the single group 
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FIG. 4. Dispersion relation in tlie valence and spin split off bands of Si calculated using double group formulation, for 14-band (solid line), 
6-band (dot-dash line) and using 6-band Weiler model (dash line). 



product bases associated with a single occurrence of a partic- 
ular single group IR. In other words, no mixing is permitted 
between states arising from different single group IRs in the 
construction of double group zone centre states. This situa- 
tion occurs only if there is no k independent inter-band spin 
orbit interaction included within the formulation. This situa- 
tion corresponds to an unperturbed Hamiltonian equal to Hq 
with the k independent inter-band spin orbit interaction se- 
lectively removed. Thus, a single group formulation with a 
limited bases set represents a system that differs from those 
described by the Hamiltonian indicated in Eq.lfTll. 

More importantly, crystals are systems that contain many 
electrons. Hybridisation is a well known many electron ef- 
fect in the formation of covalent bond in cubic semiconductor 
crystals. With consideration of spin degree of freedom, hy- 
bridised orbitals introduces additional mechanism of mixing 
within the zone centre states. Mixing due to many electron 
effects permits zone centre states to form bases of representa- 
tion of the double group, which are drawn from multiple sin- 
gle group IRs that are compatible with the double group rep- 
resentation. This mixing mechanism can not be incorporated 
directly within the one electron theory, and therefore requires 
the use of double group selection rules in the k ■ p method. 
In particular, such mixing allows the zone centre states with 
r^ and r^ symmetry to experience the effects of spin orbit 
interaction. This provides an additional mechanism that can 



modify the ordering of anti-bonding states which is absent un- 
der single group consideration. In particular, this mechanism 
provides an explanation to the observed trend that in heavier 
elemental semiconductors the Fy conduction band shifts to- 
wards the valence band. This mechanism may explain why a- 
Sn has an inverted band structure. Given that there are eight 
valence electrons per primitive cell, there are eight bonding 
and eight anti-bonding states under the consideration of hy- 
bridised orbitals. This means that each of the double group 
IRs of the Oh group occurs only once, with the positive par- 
ity IRs corresponding to valence band zone centre states, and 
the negative parity IRs corresponding to the conduction band 
zone centre states. Any other states that may exist, and are 
available for occupancy, are un-bound states above the vac- 
uum level. The effect of these states on the dispersion relation 
tends to increase the lowest zone centre conduction band ef- 
fective mass, but reduces the valence band effective masses. 
Therefore, they can not be relied on to explain the difference 
between predicted mass under single group formulation and 
that of experimentally observed values. From group theoreti- 
cal perspective, mixing in zone centre states requires the use 
of double group selection rules. The multiplicity in the de- 
composition of direct product of IRs Fg lead to difference 
in the description of A; • tt interaction between the single and 
double group formulations. 

The treatment of Zeeman interaction under the double 



34 



group formulation 'de-couples' the link between these terms 
and the angular momentum. In deriving the Hamiltonian with 
states quantised in a given orientation, the unitary transfor- 
mation required to diagonalise the Zeeman interaction differs 
from that generated by angular momentum operator Such 
transformation becomes material dependent on the k and q 
invariant for the valence band. The unitary transformation 
introduced by Luttinger only diagonalises one of the gener- 
ators, but is nevertheless a valid unitary transformation. Due 
to mixed nature of the zone centre states, the S ■ B term de- 
rived from Dirac equation has to be dealt with using group 
theoretical method. The matrix representation of S are no 
longer the Pauli matrices projected into the adapted double 
group bases but requires the general description of an opera- 
tor which transforms according to F^ IR of the Oh group. As 
generators for transformation of spin states under point group 
operation, reduced tensor elements for matrix representation 
of S is subject the same constraints as that of angular momen- 
tum operator as shown in Eq.([TO|. The decoupling of Zeeman 
interaction from total angular momentum and impact of mixed 
zone centre states on the matrix representation of the spin op- 
erator makes the theoretical prediction of the associated Lut- 
tinger invariants difficult. 

Despite the prevalence of single group formulation, group 
theoretical methods are employed to determine the number 
of invariants for a given interaction using double group selec- 
tion rules. This introduces some disparity between the number 
of independent second order interaction parameters in single 
group formulation and number of invariants obtained using 
the double group selection rules. Perturbation theory allows 
these invariants to be related to symmetric and anti-symmetric 
Lowdin interactions. The linear dependence among the in- 
variants introduced by such relations is an undesirable con- 
sequence of the the assumption of single group formulation. 
The use of perturbation theory leads to two well defined gener- 
ators for anti-symmetric Lowdin interactions. The underlying 
physics of such interactions are made clearer by the use of the 
two generators given by J and J, as apposed to J and J^ 
used in the literature, since they remove the inter-dependence 
of invariants n' and q', and highlight the contributions from 
remote states. The mixed nature of zone centre states also re- 
quires different treatment of cr ■ B term in the Hamiltonian 
in the same way as the anti-symmetric Lowdin terms. Sym- 
metry demands that the Luttinger invariants are linearly in- 
dependent from each other The double group formulation 
ensures there are the same number of independent interaction 
parameters as invariants. Therefore, the relations between the 
Luttinger invariants and second order interaction parameters 
preserves the linear independence between the invariants. In 
contrast, the relations between second order interaction pa- 
rameters and Luttinger invariants frequently lead to linear de- 
pendence among the latter under single group formulation. 

A number of models discussed in the literature have been 
compared and contrasted with the double group formulation 
of EWZ. The key features of the models, when treating 6- 
bands are summarised in table I Villi The use of selection rules 
which define the k • tt interaction, assigns each model to ei- 
ther the single or double group formulation. In addition, the 



models are differentiated according to the way zone centre 
states are classified. The double group formulation of EWZ 
on the right and the DKK model on the left represent the two 
extremes in the treatment of spin effects. 

The DKK model does not account for spin orbit interaction, 
and as such the spin degree of freedom does not arise, mean- 
ing single group selection rules apply. There are four invari- 
ants (L, M, N, and P) associated with the degenerate states 
with r^ symmetry, and four second order interaction parame- 
ters (F, G, Hi and i?2)- The model is self consistent in terms 
of group theoretical description, but does not account for any 
of the effects of spin observed in experimental measurements. 
When the bases are transformed into the adapted double group 
bases, this model gives equality among the Luttinger invari- 
ants in different blocks and the q invariant, permitted under 
double group selection rules, is absent. The relativistically 
corrected DKK model uses the same unperturbed Hamiltonian 
as the DKK model, but treats the spin orbit interactions H^^i 
and iI|o2' together with k ■ p term as perturbations. The re- 
sulting bases correspond to those of the adapted double group 
bases, and the Hamiltonian is partitioned according to the dou- 
ble group classification of zone centre states. While symmetry 
arguments show that the Luttinger invariants differ between 
blocks of the Hamiltonian, and q invariant is permitted to ex- 
ist, the evaluation of the k dependent relativistic corrections 
indicate that such differences are negligible. Thus, the rela- 
tions between the Luttinger invariants in different blocks still 
hold, and the q invariant is predicted to be approximately zero. 
Thus, the relativistically corrected DKK model fails to explain 
the finite q parameter observed experimentally, the effective 
mass of the spin orbit band, and the effective mass of the low- 
est conduction band. Under the relativistically corrected DKK 
model, the intra-band k independent spin orbit interaction is 
not relativistically small, and generally assumed to give rise 
to the finite splitting due to spin orbit interaction observed ex- 



perimentally. However, the small factor at the front of Eq.(43 1 
means a substantial number of states far away from the rele- 
vant states are required to obtain a finite result. This is perhaps 
unphysical. The most obvious reason we can call upon to ex- 
plain this is the discrepancy between the adapted double group 
bases used in the relativistically correct DKK model and that 
of the zone centre state under the double group formulation as 
discussed in section |Xl 

The Kane and extended Kane models are utilised most fre- 
quently in the implementation of single group formulation in 
the literature. The models assume double group classified 
zone centre state energies for the near set, but retain the single 
group selection rules to evaluate the k • n perturbation. It is 
also assumed that there is equality among Luttinger invariants 



in blocks associated by spin splitting (Eq.(D4f i), but the q in- 
variant is allowed to be finite. It is difficult to attribute these 
assumptions to a consistent set of group theoretical rules or 
a definitive un-perturbed Hamiltonian. The selective use of 



Eq.(D4f I and exclusion of Eq.(D4ei, both derivable under the 



relativistically corrected DKK model, suggests that a different 
classification of zone centre state energies have been used for 
the evaluation of q invariant compared with other Luttinger in- 
variants. The finite q invariant can arise by the incorporation 
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TABLE VIII. Comparison of a number of different single group formulations of the k ■ p method with the double group formulation. 



Model 


DKiP 


DKK with 
Relativistic ConectioilJS 






Extended Weiler ModeQ Ewi^ 


Selection Rules 


Single group 


Single group 






Single group Double group 


Unperturbed zonecentre energy 


E(r+) 


E{rt) 






E(r+):_E(r+) E{Tt);E(T+) 


Predicted zone centre energy 


s(r+) 


B(r+);E(r+) 






E(T+)-E{T+) E(Tt)-.E(T+) 


Invariants 


L, M, N, P 


7l,72,73,K,g = 
7r = 72, 7^" = -73, 

K - K g- 

, li" = 7i,«'' = 4k + go 






f 71, 72,73, K,« f 71,72,73, K,g 

{ 7r,7r,'«''° { 7r,7r,'*"° 
1 7r, «^°° 1 7r,«°° 


Second order 
Interaction Parameters 


F, G, Hi , H2 


F,G,Hi,H2 






^t r+,r+ 1 r+,r+i 

7^8 ;i ^rg ;2 
P+' ^p+ p+' ^p+ p+ 

1 7 i 8 •' 7 ^ 8 '^ 7 
2^8 

r+' r+,r+ 


b 


■ 


r+ r+' r+ r+ ' 
' 8 ''s ^ 8 ■' 8 

Fg ;1 Fg !2 Pg ;3 

V+ p+' '»p+ p+' V+ r+ 

'8 ''S ^ 8 •' 8 ' 8 ''S 
^1-7 ,rg ;l Pg ,2 

.r-v rg 


Bases 


Single group 


Adapted double group 


bases 


or sing 


e group product bases Double group bases 



'^ The model was originally described in Ref. 17 and extended in this manuscript. 
'" This parameter is only present if remote state with V7 symmetry is included. 



of States with Fg F^r (FJ ) symmetry into the near set in the 
extended Kane model. If the same energy scaling is applied 
to the evaluation of other Luttinger invariants, it is clear that 
relations in Eq.( |D4f| ) are also invalid. Therefore this calls for 
double group classification of all zone centre states, as in the 
Weiler model. 

The Weiler model classifies all zone centre states accord- 
ing to the double group, but retains the single group selection 
rules when evaluating the k • -k perturbation. While there is 
no mixing of states between different single group IR, as in all 
single group formulations, this model constitutes an improve- 
ment over the relativistically corrected DKK model. This is 
principally due to the double group classification of zone cen- 
tre states in both the near and remote set which ensures that 
Luttinger invariants differ between blocks of the Hamiltonian. 
In the extended Weiler model, discussed in this manuscript, 
there are ten independent second order interaction parameters 
and ten invariants describing a 6-band model. Ten indepen- 
dent second order interaction parameters exist only if Fg re- 
mote states derived from F^ IR are included in the near set in- 
teractions. In the original Weiler model, the Fg^ remote states 
derived from Fg^ IR are excluded. Parametric fitting of the 
full set of second order interaction parameter will also lead 
to un-physical values for some of the interaction parameters 



7r8,3 

'r+,r;J 



> 



7^8 '2 



). The inclusion of additional conduc- 



tion band zone centre states will continue to have difficulty in 
explaining the conduction band effective mass in this model. 

The double group formulation of EWZ follows naturally 
the prescription of Lowdin method in partition of states into 
near and remote sets. Compared with the Weiler model, mix- 
ing of double group bases derivable from the compatible sin- 
gle group IRs is permitted in addition to the double group 
classification of the zone centre states. From a one elec- 
tron theory perspective, this takes into account the k inde- 



pendent inter-band spin orbit interactions among the remote 
set of states. The double group selection rules also enable 
the semi -empirical method to account for the effect of mixing 
due to formation of hybridised orbitals with consideration of 
spin. This model has the same number of independent sec- 
ond order interaction parameters as invariants in each block. 
Furthermore, the same set of second order interaction param- 
eters are compatible with effective masses observed experi- 
mentally in the lowest zone centre conduction band, valence 
band, and spin split-off band. Compared with the extended 
Weiler model, the un-physical values of the second order pa- 
rameter can be explained through mixing, which occurs in 
all states with Fg symmetry. Therefore, the number of in- 
dependent second order interaction parameters in each block 
is already the same as number of invariants without the need 
to specifically include remote state derived from Fj^ symme- 
try. All the other positive features of the Weiler model are 
retained, enabling good description of the electron dispersion 
near the zone centre. 

The Mj quantum number and its correlation to the z com- 
ponent of angular momentum have been used extensively in 
the explanation of many physical observations such as circu- 
lar dichroism. However, the historically chosen unitary trans- 
formation, which diagonalises the k independent intra-band 
spin orbit interaction, leads to an order of basis in the valence 

band zone centre states which is not the same as Dt IR of 

2 
the 0(3) group. Hence, one can not make a correlation be- 
tween the Mj number of the historically chosen basis and the 
z component of angular momentum. The difference between 
the historically chosen basis and those bases with the same 

order as D\ IR of the 0(3) group, can easily be accounted 

2 
for by a similarity transformation. As a result of this trans- 
formation, discussed in section [rV] the form of the inter-band 
Hamiltonian in a 6-band model differs. Also there are some 
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sign differences between the Luttinger parameters, most no- 
tably in the change in sign of 72 invariant in the valence band. 
Under the double group formulation, the matrix representation 
of the angular momentum operator consists of two linearly in- 
dependent matrices, as shown in Eq.(|9]l. The explanation of 
phenomena, such as circular dichroism, can still be made pro- 
vided appropriate momentum matrix elements are evaluated 
between relevant states using group theoretical techniques. A 
prime example of this, is their applicability in deriving the se- 
lection rules on polarisation dependence of optical transitions 
between valence and conduction band zone centre states, as 
discussed in sectionlVl 

Experimentally determined Luttinger parameters are de- 
pendent on the model used to describe the interactions be- 
tween the electro ns and the external perturbation such as stress 
or magnetic fielcP^El Specifically, they are dependent on 
the model derived under the single group formulation, which 
utilises the relations given in Eq.(D4f 1. The experimentally 



obtained Luttinger parameters are then generally not in agree- 
ment with the relations derived between the second order in- 
teraction parameters and Luttinger invariant under the double 
group formulation. The inter-band Luttinger parameters have 
not been determined independently either There is a need to 
re-analyse the raw data and obtain new Luttingers parameter 
using models which follow double group selection rules. The 
difficulty in obtaining inter-band Luttinger parameters is then 
partially overcome, if there are measurements of Luttinger pa- 
rameters describing relevant intra-band interactions. Such pa- 
rameters, have a strong influence upon the dispersion away 
from the zone centre. With the advent of modern photoemis- 
sion instrument92323 a case can be made to re-study these old 
material systems in order to obtain appropriate Luttinger pa- 
rameters. For the same reason, the parameters shown in Table 
|VII| and dispersion relation obtained under the double group 
formulation, shown in Figure 3]4 may not be reliable due to 
its dependence on Luttinger parameters extracted using single 
group formulation. 



states belonging to a particular IR with dimension larger than 
one are degenerate, linear combinations of such states are also 
energy eigenstates, and form basis of the IR. One is free to 
choose such basis, but the basis with the same ordering within 
in the IR must be used for calculation of all the relevant inter- 
action matrices. This applies to first and higher order terms in 
perturbation theory. 

It is customary to generate representation matrices by con- 
sidering the point group as a subgroup of 0(3). Using the 
total angular momentum operator J as generators of the in- 
finitesimal rotation in the 0(3) group, representation matri- 
ces for proper rotations can be constructed, and effects of im- 
proper rotations in the Oh or Td groups can be introduced 
afterwards^^. For the Oh point group, the effect of improper 
rotations differentiate representations that have different par- 
ity. These representation matrices fix the order of states within 
an IR. This aspect is essential in deriving the necessary simi- 
larity transformation required in combining the basis of Fg in 
the Oh group, into basis of Fg in T4 group under compatibil- 
ity relations. If a different basis is used, then the underlying 
unitary transformation must be considered. It is preferable to 
work with system with higher symmetry (Oh group in this 
case) and consider the effect of lower symmetry {T^ group) 
via the sub group relationship and compatibility relations. For 
a Hamiltonian constructed using the method of invariant^, the 
form of the matrix representation of components of the total 
angular momentum operator with respect to the J — "^ va- 
lence band bases, implies an order of the states same as those 
used in EWZ^^. 

It is also important to consider the time reversal properties 
of basis when considering intra-band interactions. This is of 
particular importance when obtaining first order, intra-band 
interactions, allowed under Td group symmetry. To facilitate 
this, we use the relation stated by Kostep^and used in EWZ, 
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Appendix A: Basis functions used in calculation of interaction 
matrice 

The solution of un-perturbed Schrodinger equation (zone 
centre states) form the basis of irreducible representations of 
the symmetry group of the crystal^^ . The first order interaction 
matrices between states belonging to irreducible representa- 
tions may be calculated using the extension of Wigner-Eckart 
theorem for the point group^. This reduces the interaction 
matrices to a sum of products between a scaling parameter 
(reduced tensor element), and an angular dependent part. As 



where T is the time reversal operator, \J, Mj) is a basis of IR 
of dimension 2 J + 1 and Mj = —J, —J + 1, • • • , J index 
the basis within the IR. We ensure the double group bases 
follow this convention, but retain the real harmonic format 
generally used for single group bases. However, a necessary 
scaling factor of i is included in the even single group bases 
(expressed as real spherical harmonics) to ensure compliance 



to ordering of basis prescribed by Eq,(Al 1. This also ensure 



the reduced tensor element of momentum operator is real. 

Table |IX] lists the basis functions of representation of the 
Oh group used in the calculation of interaction matrices in ap- 
pendix IB] The single group bases are expressed in terms of 
real harmonics whereas the double group bases are expressed 
in terms of real harmonics and spinor states {a, /?}. The dou- 
ble group bases are compUant to the time reversal condition 



stated inEq.lAl 1. 
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TABLE IX. The single group bases are expressed in terms of real harmonics, whereas the double group bases are expressed in terms of real 
harmonics and spinor states {a, /?}. The double group bases are compliant to the time reversal condition stated in Eq.l|Al|. 



IR 



r+ 

^ 3 



r^^ 



basis functions 



IR 



basis functions 



r+ 

^ 7 



1 = 0) 
1 = 2) 

1 = 2) 

1 = 3) 

= 0) 

1 = 2) 

1 = 2) 



6z — r ^ 

V3{x'-y') 

yz 

zx 

xy 

x{3y^ + 3z^ - 2x^) 
yiZz^ + 3^2 - 2j/2) 
2(3a;2 + 3j/^ - 2z^) 

ia 

i/3 

xya + (yz + izx)f3, 

—xyj3 + [yz — izx)a 

— (yz — izx)a — 2xyP 
VS{yz — izx)l3 
—y/3{yz + izx)a 
(yz + izx)j3 — 2xya 



/c 2 2 2 

yz(bx — y — z 



{1 = 4:): { zx{6y^ 



xy{6z — X — y 



il = 2):{ 



i{x^-y^)P 
^Jr^-3z^)a 

-i(x'^ -y^)a} 



r4 



rr 



Tt 



To- 



1 = 3): xyz 

V3xyz{ 



1 = 1) 




2 2\ 

X -y ) 



xyz{r^ — 3z^) 



x{y'-z'), 



-' y') 



viz' 

z{x^ 

ixyza, 

ixyzP 

— za — [x + iy)P, 

zP ~ [x — iy)a 

~{x + iy)a, 

I [22a ~{x + iy)l3] , 
I [22/3 + {x- iy)a] , 
I (x - iy)f3 



Appendix B: First order interaction matrices 

The perturbation operator, -k = p + ,, \ 2 \S x Wir) ) 

(or p under the single group) transforms according to r4 IR of 
the Oh group, and is time reversal odd. A general irreducible 
tensor operator is time reversal even. Thus, the equivalent op- 
erator of TT constructed from irreducible tensor operator must 
be scaled by a factor of "i". The angular dependence of -k with 
respect to given bases, may be evaluated using the constructed 
equivalent operator, Wigner-Eckart theorem. Therefore, 

7^—1 7^+1 

^ 1 1 

^x — 7= 

rp—\ I rri+l 



where Tl are irreducible spherical tensor operators of rank 
one. Given the set of zone centre states associated with the 
a" , a" IR of the symmetry group, the first order interaction 
may be written in matrix form as. 






mo 

h 

mo 



Aie{a;,y,z} 



'^fk^l'/'? 



(B2b) 



Where selection rules indicate that there are multiple linearly 
independent angular parts, they can be generated by enumer- 
ation of basis of the IRs. 

With the single group bases given in appendixJAJ the matrix 
representation of the k • p interaction can be written as. 



i^P+ p- 



— e 



mo''r+,R 



1 r\irf~ ryirti K 



K, 



p+ p- 
^ 5 '^ 4 






mo^r+.r^" I ^^ k: 




ma ^r^ ,r+ ykx ky fc. 



r-,r+ - mo^r-,r3 



Kj.+ p- 




(B3) 
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With the double group bases given in appendix IA] the matrix representation of the fc • tt interaction can be written as. 



/C 






r±.rT - mo5r±,rgT I , 






^r±,rj - mo^rJjgT 



-Art. 



JC 






= /C 



J^ a iJ^ a 



V V3A;+y 



JC 



^ a 



-*- a i-*- a 



mo 



A a ?i a 



^r±,rT - W^r±T,T 



/fc, fc+\ 

-fc, fc_ 

fc+ fc, 

Vfc_ fc,/ 



/ k+ 2fc, \ 

-V3fc+ 

ysfc- 

V 2fc, -fc_ / 

/ 4fc, ysfc- 
\/3fc+ 3fc- 
mo^^8^i8 3fc+ 



;^2 

?p± pT 



(B4) 



fc+ \ 

/3fc. 
V fc_ V3fc+ -4fc,y 



(B5) 



r 



To obtain the individual component matrix M^^^ ^„ , one takes 
the corresponding Ka^^av and set fc^ = 1 and fc^ = OVz^ ^ ji. 
The first order interaction with respect to the double group 
bases may also be obtained by transforming the first order in- 
teraction matrix with respect to the product bases into adapted 
double group bases, using the transformation matrices given 
in the following section. This is illustrated in appendix IE] In 
applying Winger-Eckart theorem to obtain these interaction 
matrices, the first order interaction parameters ^'s are real, but 
their signs are determined by the radial part of the wave func- 
tions. 



Appendix C: Transformation matrices from single group to 
adapted double group bases 



The Oh point group character table indicates that Fg 



r7 = r; 



and r^ 



r^ 



4 — -^ 8 ^ -^ 6 ' """ ^ 6 ^ -^ 5 

the adapted double group basis for Fg 



- F+ 



Fy . Therefore, 
F^ and F+ F+, 
may be obtained from a direct product of spinor states with 
single group basis of F4 and Fjj", followed by an appropriate 
unitary transformation. All bases transform as row vectors. 
The unitary transfomiatiorP^may be written as. 



(f^^ , ^l^ , ^l^ , Vl^" , ^l^ , ^l^ ) = (c^l^a, 0: 



^ p+ p+ p+ p+ p+ p+ 



(^F- 



{j>i' a, 



2 '^' 



,4^/3 



1 

_ i 







V 







2 

V6 
1 

V6 

_ i 

V6 





1 

V6 

_ i 

Ve 






2 

Ve 



I 76 



,0?/3 



\/6 





2 

Ve 








1 

_ i 





1 

\/2 



V2 








(</>[" a, '^'2" "> '/'a' "' '^i'" /3' 4'1'' /3, 03' 13 j U^ 








1 

_ i 





V6 








1 

\/3 

1 

_ i 
'\/3 





h\ 



76 73 

1 1 

\/6 \/3 



\/3 





73 

i 

73 





1 

73 



73 

73 




y 



73 



(Cla) 



(Clb) 



The two matrices are related by 



t/'4 =[/i 



L/ 




(C2) 



where L/ is defined by Eq.( 12 1 
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Appendix D: Construction of the Lowdin term and relations of Luttinger invariants to interaction parameters 



Without consideration of spin, the DKK hamiltonian (i?L5tud,;n) with respect to the V^ basis of {yz, zx, xy} in the valence 
band, are obtained from first order interaction matrices, and Eq.(17c|17d|i as. 



Lkl + M{kl + kl) Nk^ky+iP[k^,ky] Nk^k,-iP[k,,k^i 



Hdkk — I Nkxky 



iP%M 



L = F + 2G 
M = Hi+ H2 
N = F + Hi-G-H2 

P=\{F-Hi-G + H2 



Lkl + M{ki + kl) Nkyk, + Pi[ky, k. 



1 V Kij f^z 



iP[ky,k,] Lkl + M{kl+ki) 



(Dla) 

(Dlb) 
(Die) 
(Did) 

(Die) 



where the second order parameter F,G,Hi,H2 are as defined 
by Kane'^. These second order parameters may be related to 
first order interaction parameters in appendix[B]by, 



F 



G = - 



Hi 



H2 = - 



6cr = 

27770 


h^ 


^r+,r- 




ml El^ - Ef 








2 


27770 


h^ 


%.r- 




2ml Sp+ - ^p- 








2 


67r = 

27770 


h^ 


^r+,r- 




ml E- - Ef 
'-5 J- 4 








2 


^%.- 


te 


%^r- 





27770 



^oE^^ 



E. 



f 



(D2a) 



(D2b) 



(D2c) 



(D2d) 



where the sign of each term in the summation is determined 
by the relative energetic position of the remote states and near 
states. For diamond. Si, and Ge, all are negative. The parame- 
ters, a, 6, IT and (p are dimensionless versions of the second or- 
der interaction parameters for the valence/spin split-off bands 
defined along the lines of Refi50landl78l 

The 6 band Hamiltonian with respect to the adapted double 
group basis can e derived from obtained from 



variants and second order interaction parameters as, 

27770 (i + 2A/) 27710 {F + 2G + 2Hi + 2H2) 

7l = n; ; 1 = -^1^ : 1 
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73 



^72 3 h^ 3 

2(7 + 47r + 4(5 + 40 - 1 

2mo (i - M) 27770 {F + 2G-H1- H2) 



h^ 6 




h^ 
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— (T + TT — 


26 ^(f) 








27770 A^ 


27770 {F - 


-G + Hi- 


-H2) 
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a + TT — 5 


-4> 
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50 


27770 {F-G- 


H1+H2) 





;72 3 6 fi2 


6 
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6 
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ir = 72 


73 = -73, t + y " 


, .90 

= '^ + y 


ss 

7i =71 


s 5o . z' , 5o\ 
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These results have been obtained previously in the literature^. 
Within the framework of single group approximation, there 
are four interaction parameters (F, G, iJi, and H2) and four 
invariants (L, Af , N, and P) permitted by single group se- 
lection rules. Referring to the adapted double group bases, 
the single group selection rules only permit four invariants 
(71 1 72, 73, and k). Note that the k invariant in the spin-split 
off band differs from that of the valence band. Under this par- 
ticular approach, the effect of inter-band k independent spin 



H-^-P = H^ + U'a {t2y.2® HdkkW^ +gofiBS-B 

(D3) 
where the U^^ is given in appendix O This is formally the 
DKK Hamiltonian, with k independent relativistic correction 
only. The effect of Hso3 = QoIJ-bS ■ B is also included by 
projecting the matrix representation of S into the appropriate 
spaces. It returns the form of Hamiltonian given by Eq.( 111, [Bland Eq.( 17c 1, the symmetric Lowdin term has the form of 



orbit interaction, described at end of section III and beginning 
of section |IX] could not have been included. 

In the double group formulation of EWZ, the second order 
interaction parameters (^'s are defined by Eq.(ll-13) of EWZ. 
Using the first order interaction matrices obtained in appendix 



(D4a) 
(D4b) 
(D4c) 

(D4d) 

(D4e) 

(D4f) 



with the Luttinger parameters related to the single group in- Eq.([TT]i with the Luttinger invariants defined by. 
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71 = 2(11 j,+ + 2di , + dlL + sdrl 

-^8 '-^8 -^S'-^S -^8 '-^8 -^8 '-^8 









73 — Ct^+ r+ + Ct^+ r+ ^Cr+ r+ 









Sp+ p+ 

^ 8 i'- 7 

Sp+ p+ 

^ 8 i'- 7 



7i 



A^7 

''r+,r+ 



4Cp+p+ 



7i^(r7) = c;^^-+4C 



r+ 

^ 8 

r-,r; 



-1 
- 1 



4c;rp+ 

i 8 -J- 7 

sCptl 



■4Cpr;^+ 

'^ 8 i'^ 8 



(D5a) 
(D5b) 
(D5c) 
(D5d) 
(D5e) 

(D5f) 

(D5g) 



r 



It is clear from Eq.( D4b|D5b i that the sign of 72 is determined 

by relative strength of <t versus n under single group formula- 

p^ p^ 

tion, and CI „+ versus T + p+ under the double group formu- 

^ 8 ■ 8 ^ 8 ■ 8 

lation. Given that the definitions of these parameters contain 
the energy differences between the corresponding zone centre 
states, the sign of 72 is mostly determined by the symmetry 
of lowest zone centre conduction band states with the possi- 
ble exception of Si. Thus for most materials except diamond, 
and possibly Si, 72 < for the form of Hamiltonian given 
in Eq.dl 1. In Si, the contribution from other parameters in 



Eq.( D4b|D5b i may become important in determining the fi- 
nal sign of 72. While intra-band absorption data suggest a 
positive 72'*^, cyclotron resonance data under uniaxial stress'*" 
suggests a negative 72. 

Group theory based methodologies are empolyed to obtain 
the first order interaction matrices in appendix |B] but can not 
specify the sign of ^ parameters. The second order interaction 



parameters ( are all positive for the main valence band block. 



except the mixed term of <^. 



rr,3 



The sign of Luttinger pa- 



^P+ p+- 

'- 8 ''- 8 

rameters can, therefore, be inferred from the relations shown 
in Eq.(D5|23 I obtained using perturbation theory. 



Appendix E: Relationsip between first order parameters 

The first order interaction between states belonging to any 
single group IR have been listed in appendix IB] Using appro- 
priate transformation matrices between the product bases and 
adapted double group bases (see appendix |C]l, the interaction 
between states of F^ and Tj referring to adapted double bases 

isgivenby^p^^Pg:P^,g,Pg = U^' (lL2x2«) A'p^^pjt/^^. Ifwe 
consider the interaction between F^ and F^ , we have. 



'^r+er+(r+),r, 



O^eC^i) 



= U^^ 



+ t 



1-2x2 ' 



mo 



'F+t: 




'- 5 ■'- i 
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V2 
2fc. 
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fc_ 
Ve 

2fej 

V6 
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(Ela) 



(Elb) 



(Elc) 



where fp+ p- = P^ ■ Comparing the reduced tenor elements of each block, we obtain: 



pp 
p 



e 



rt(p)'^i(p') 



"^^r-rr7 



-^^r+(r+).r-(r-) 

= \/6Cr+(r+),r-(r-) 

= \/6^p+fp+i.p-fP7i = ^Cp-fr7 



r8"(r4"):r+(r+) 
^^r3-(r-),r+(r+) 



'r^(r+),rg-(r- 



r6"(r4),r+(r+) 



(Eld) 
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The interaction matrices refer to the adapted double group ba- 
sis given in appendix |AJ with the symmetry of the IRs deter- 
mined from the character table. It is worth noting that the rela- 
tion between signs of first order parameters are fixed between 
different blocks, even though the symmetry does not permit 
the determination of the actual signs. Particularly, the sign of 
£_ in /Cp+,p+, p-fp-) is negative relative to /Cp+,p+, p_ 



Similar relations in Eq.(45 i may be obtained under the sin- 
gle group formulation by transforming the appropriate inter- 
action matrices that refer to single group product bases. 



Appendix F: Group theoretical background behind the unitary 
transformation of Eq.(12) 

When applying group theoretical method to problems such 
as crystals with diamond and zincblende lattices, where a 
group and subgroup (Tii C O^) relation exists, it is advan- 
tageous to study the higher symmetry Oh group, and examine 
the subgroup under compatibility relations. This approach has 
been taken by EWZ, and has been used successfully to derive 
all necessary generators for the Td group, including those that 
are spatially odd, such as the linear k term, from results of 
the Oh group. It is not generally possible to take this process 
in reverse and infer the properties of parent group from the 
structure of the subgroup. 

Under the method of invarianP, the generators for irre- 
ducible perturbations are obtained from known bases of ap- 
propriate IRs associated with the zone centre states. The form 
of these generators, as well as the form of the representation 
matrices, are dependent on the choice of bases. The point 
groups, such as Oh and Td, are subgroups of O (3) . As a result, 
the IR of the 0(3) group also form representations of the point 
group which are generally reducible. For the specific cases of 
Dj and Df IRs of the 0(3) group, they also form irreducible 

2 2 

representation Fg and Fg IR of the Oh group, respectively. In 



the literature, the AIj number used to label basis of IR of the 
point group is frequently equated to z components of angular 
momentum. It must be recognised that in point group, there 
is no corresponding generators of infinitesimal rotations, as 
in the continuos Lie group such as 0(3) or 5*0(3). There- 
fore, no link can be established between the Af, number and 
z component of angular momentum in the same way. If one 
is to correlate the Mj number as scaled z component of angu- 
lar momentum, then the order of bases must be derived from 
the appropriate representation of the 0(3) group. In the 0(3) 
group, the infinitesimal generator of rotation about the z axes 
is scaled from z component of the total angular momentum 
operator The order of bases within IR Dj^ has the rrij num- 
bers ranging from 2j + 1 to — (2j + 1), which labels the row of 
the representation matrices. The corresponding order of bases 
in the IRs of the point group are then given by the correspond- 
ing \J,Mj) numbers that correlate to the \j,mj) number in 
the Df representation of the 0(3) group with j = ^ for Tq, 
and j = § for F± in the Oh gi'oup. For the Df IR of the 0(3) 
group, the representation matrices are the same for proper ro- 
tations but differ in sign for improper rotations. Within the 
Oh group, such differences in the representation matrices are 
reflected in the isomorphism between the product of IR rep- 
resentations F7 



F- and IR r+ (or F+ = F^ 



F-). The 



corresponding representation matrices are hnked by, 
D^t (5) = i^r- (^) ^ ^r- (^) ^ ^rr (^) . ^r" ^^^ ^p^^ 

since T^ is one dimensional and £''"1 (g) is just a number 
The representation matrices and corresponding bases func- 
tions derived by Onodera et al. follows the relations shown 



inEq.(|FTi. 

Similar relations between IRs may be constructed from 
other one-dimensional representations of the Oh group. One 
such example is 



Ff = Fo 



r- {t, 3} or {j, i} e [{1,2}, {3, 3}, {4, 5}, {6, 7}, {8, 8}] 



D^t'(g) = i^r- (^) ^ ^r; ^^^ ^ ^r- (^) . ^r; ^^^ 



(F2) 
(F3) 



If one makes the same choice of basis, and representation ma- 
trices for Fg^ as in D^, then there will be some sign differ- 

ences between the corresponding representation matrices of 



D^s (g) (from Eq.lFl 1) and D^s'{g) (from Eq.(F3 1) depend 



ing on the value of D 1 (g) and D 2 ((^). Both clearly form 
representations of the F^ IR, and must be linked by a similar- 
ity transformation U, which relate the basis functions. 



D^s'^g)^U^D^S^g)U. 



It can be shown that U is defined by Eq.(|T2|l. 
In the case of a general matrix of the form 

p — ix<n I p „ U" \ 



(F4) 



(F5) 



I 

where the operator P transforms according to F^m IR of the 
symmetry group, its transformation properties are given by the 
product representation F* „ ® Fq,™ (g) F^i., which is generally 
reducible and decomposable into IRs. The matrix represen- 
tation P is represented by a number of linearly independent 
matrices and associated invariants depending on the decom- 
position of the product representations. Specifically, the num- 
ber of linearly independent matrices is given by the multiplic- 
ity factor of Ta^ in the decomposition of F^„ (g) Va^^^ ■ To 
find these linearly independent matrices, one must utilise ir- 
reducible tensor operators transforming according to Fa^ IR 
and bases of IR F^u and T^^ . This procedure has been used 
in EWZ utilising the bases given in appendices iBlto calculate 
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the first order fc • tt interaction matrices. 

Since F^* ® Fj^ — V\ and V^ is one dimensional, we can 
write, 

F;„ ® F„» ® F„. = (F^* ® F;„) ® r„™ ® (F^ ® F(P^) 



and the corresponding relation in matrix representation may 
be written as 



D^''^- (g) ® D^'- [g) (g> D^^^' (g) = (^D^^' (g) ® D^'^^ (5)) (g> D^^^ (g) (g> (^D^^ (g) ® D^^^ (g) 

= {d^^' (3) ® D^'"- (5)) ® [d^- (g) ® D^'- (.g)) ® D^-" (g) 



(F7a) 
(F7b) 



r 



Historically, the point group bases are not well known, and 
calculations are generally made utilising the fact that the bases 

of of and DJ IRs of the 0(3) group also from represen- 

2 2 

tation of the Fg and the F^ IRs of t he O h group. Trebin 



et al.t^has utilised the relation in Eq.(F7ai to show that the 
generators between states of any of the {F^, Fy', F^} triple, 
are the same as the corresponding generators between states 
of the {r([,Fg ,Fg"} triple. Thus, the full 8-band Hamil- 
tonian for the zone centre states transforming according to 



should be considered with introduction of appropriate fac- 
tor of "i" (general irreducible tensor operator is time rever- 
sal even, angular momentum is timer reversal odd). For 
example, the generator transforming according to F^™, be- 
tween states of the Fg and F^ IR, is the same as the gen- 



erator of Fo 



between states of the Fo and Fo IR 



scaled by a factor of "i". Thus, all the generators may be 
obtained from the {F^i^ , Fg , Fg^ } triple. If we consider the 
generator for the tv operator transforming according to FJ 



{Ff ,r+ r+} triple can be constructed from generator ob- ^^^^^^^ ^^^^^^ °f ^s and Fg symmetry, then they are the 
tained from the {F+, Fg", Fg } triple. This leads to a form of ^^"^^ ^' ^^^ °f generator for F^ - - - -+ — — 



Hamiltonian similar to those shown in Eq.([T5]l. Such a pro- 
cess does not require any knowledge of point group bases of 
the {F^,F^,F^} triple. The methodology is clearly valid, 
but should be used with the following considerations, 

1. the vector space of F^ IR is orthogonal to that of Fg^ 
IR, and the bases of Fr IR are not the bases of Ft IR. 
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the method is not applicable if either F^k or Fq,i. does 
not belong to the {Ff , Fy^ , F^} triple. 



the representation matrices generated by D^^'{g) = 
D^s (g) (g) Z^rj (g) = £)^s (g) . D^2 (g) constitute one 
of the many equivalent representations. It is equ ivalent 
to the one commonly used for symmetry analysis^^ElEll 
by a similarity transformation shown in Eq.(|F4|l. 



4. correlation of the Mj number in the D^^'{g) repre- 
sentation with the z component of angular momentum 
is not appropriate, since it is not compliant to relation 
A unitary transformation of the basis described 



Eq.lFl 



by the inverse of Eq.([T2]i is required before such corre- 
lation can be made. 

The second point leads to the scarcity of information in the 
literature on the derivation of the two linearly independent 
first order interaction matrices, between states of Fg IR. This 
deficiency is overcome by considering the relation shown in 



Eq.(F7bi. The generators with F„^. symmetry between states 



associated with IRs in the two different triples is the same 
as that with F^ ® F" symmetry between states of com- 
patible IR in the one of the triples. The compatible IR in 
the opposing triple is obtained by direct product of the orig- 
inal with V2 ■ The time reversal properties of the operator 



— F^ between 
states of Fg^ and Fg^ . These generators have been obtained as 
matrix representations of i{Jy, Jz], i{Jz, Jx}, i{Jx: Jy} and 
{J,, ( J2 - Jl)}, {Jy, (J2 _ j2)|^ I j^^ ( j2 _ j2)| ^ith re- 

spect to the DJ basis. These generators refer to the basis set 

2 

for which D^>i '{g) are the representation matrices of the F^ 
IR and the Mj number does not correspond to the z compo- 
nent of angular momentum. The second set of generators has 
not been use d in the co nstruction of multiband Hamiltonian in 
the literatur^nHHHi. 

The unitary transformation T introduced to the basis of Fg^ 
in EWZ, ensures that the compatible Fg^ and F^ IRs have 
the same representation matrices as F^ for the elements of 



the Td group. This corresponds to Fg — Fg 



In 



essence, this is the reverse of the process that infers repre- 
sentation matrices of F^ IR from that of F^ IR. Applying the 
inverse unitary transformation to the F^ bases, yields the re- 



sult obtained by Chuang as shown in section IV Thus, one 
would expect the transformation matrix inE q.([l3]l to be given 



by M = U^i If^i as defined in Eq.(|C2[). There is, how- 



ever, difference in the sign of the basis for the F^ and Fg IRs 



used, which accounts for the difference between Eq.( 13 1 and 



Eq.(C2i 



Appendix G: Unitary transformation induced by rotation of 
coordinate system 



The diagonalization of the specific magnetic interaction in 
an arbitrary orientation is always possible for any zone centre 
states belonging to the F^. F^, and F^ IRs. The process is 
relatively simple for states with Fg and Fy symmetry. This is 
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due to the one linearly independent generator required under 
double group selection rules for operator or irreducible per- 
turbation that transforms according to T^ IR of the O^ group 
for states belong to one of these IRs. On the other hand, two 
linearly independent matrices are required when we consider 
matrix representation of Zeeman interaction between states 
with r^ symmetry. Thus, the two linearly independent ma- 
trices of the third component of angular momentum operator 
can not be simultaneously diagonalised in general (they can be 
simultaneously diagonalised only when referring to the {xyz} 
coordinate system) for zone centre states belonging to the Fg 
IRs. This is one of the characteristics of the double group ba- 
sis that can not be represented by a direct product between 
single group basis from one IR and spinor states. If one really 
wishes to diagonalise the matrix representation of the Zeeman 



interaction, in the direction of the field, the unitary transfor- 
mation should be those that diagonalise the total Zeeman in- 
teraction kJ^^^ + qJ^^^- In this case, a single diagonal ma- 
trix would be obtained. This is in contrast to the suggestion of 
Luttinger** where the transformation described by Eq.(47) of 
Ref 8 was assumed to have re-quantised the zone centre states 
along the 3 axis. In fact, it was the J'^^^ generator that was di- 
agonalised but the Jg^'^ generator remain non-diagonal. Here 
we will continue the treatment by Luttinger, and use J' as the 
generator for the induced unitary transformation on the basis 
of r^ IR. This leaves J^^^^ diagonal, but J^^ non-diagonal 
and orientation dependent. 

The unitary transformation which diagonalised J'-^^^ in co- 
ordinate system {123} are defined by the appropriate matrix 

representation of the angular momentum operator scaled by 

1 

h' 



U^ = U1U2U3 = cxp[-iGla] cxp[-iG^,/3] exp[-iG'^„7] = exp[-iG • 6] 

I 



(Gl) 



where G^ has the form of Pauli matrices scaled by i for Fg , 
and F7 IR, J^ for Fg^ IR, and J^ for F^ IR, a, /3 and 7 are 
Euler angles which define the rotation 6 that takes the {xyz} 
coordinate system into {123} coordinate system. J^ and JT^ 
are those used in by Eq.(j9]l. The similarity transform for the 
Fg transformation are related by. 



For magnetic field aligned in the [111] direction (the 
{1,2,3} axes are in the [Tl0],[112] and [111] directions), 
a = — 7r/4, /3 = — arctan(\/2) and 7 = 0. The unitary 
transformation are given by. 



U^8 (6>) = rt^r, 


' {e)T where (G2) 


/ 1 o\ 


T = 


0-1 
-10 


(G3) 


\ U i u u / 


U^^rnni])-^ ( (« + 6*)cos(|) (m-6)sin(f) 
^ ^^ ^' 2 \-{at + b)s[n{l) (a-6z)cos(|) 


C/r.-([lll]) = i 


/ 4cos(§)'^ (ai + 6) 
^/3(6i + a) (sin(f ) + sin(M)) 
y3(cos(f)-cos(f)) {hi- a) 

\ -4sin(f)^ (&i + a) 






-^ 


(cos(f)-cos(M))(ai4 



^/3 



sm 



'30 



)) (M-a) 



(cos(f)+3cos(M)) (6i + a) 

(sin(f)-3sin(^)) (&i - a) i 

^(cos(f)-cos(^)) (ai-6) 



b) 



'39 
2 
v3 



4sin(f) {hi -a) 



i(f ) - 3 sin(^)) {hi + a) i -73 (cos(|) - cos( 



3e^ 



^/3 



-4cos(|) 



/ 



(G4a) 



<-^ 



- (cos(f ) + 3 cos(^)) (6i - a) -^3 (sin(f ) + sin(^)) {hi - a) 



U^s ([111]) 



— cos 



M- 



- (cos(f ) + 3 cos(^)) (M - a) 
-^/3(M + a) (sin(f) + sin(f )) 
^))(ai + 5) 
^))(ai-6) 

M))(6i-. 
{hi + a) 



\/3 (cof 
\ — (sin(l) — 3 sin 

-73 (cos(f ) 



-4sin( 



cos 

\3 



-73 (sin(§) 
-4cos( 
-4sin( 

/3 (cos(f ) - 

- (-n(i) - 



sm 
n3 



■30 
. 2 



(ai - b) 
)) (ai + b) 



(G4b) 



2y (ai — 6) 
Ifibi-a) 
cos(^)) (6i 



<-^ 



4cos(|) {ai + h) 



3sin(M)) (ai + h) \ 
-73 (cos(f)-cos(^)) {ai-b) 
-73 (sin(§)+sin(M)) (bi-a) 



(G4c) 



-73 (ai - b) (sin(l) + sin(M)) (cos(|) + 3 cos(f )) {hi + a) J 
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where a = \/2 + y^, b = ^2-^2, 6* = -/?, and tan(6') = ^2. 

For magnetic field aligned in the [110] direction (the {1,2,3} axes are in the [001] , [110] and [110] directions), a = — 7r/4, j3 = 
— 7r/2 and 7 = — 7r/2. The unitary transformation are given by: 



t/^6.^([110]) 



C/^s ([110]) 



C/r« ([110]) 



1 



— 6 ai — b 



2(1 + i) \bi-a a-bij 

I — ai — b — \/3(ai + 6) — vSai — ^/ib — ai — 6 \ 

1 — i — v3(— 6 + ai) 6 — ai ai — 6 v3(&ii + ai) 

8 \/3a — ^/?>b\ b\ — a bi — a \/3a — y/3bi 

\ bi + a — v3&i — v3a v36i + \/3a —bi — a y 

/ 6i — a VSbi — v3a vSbi — v3a fei — a \ 



1 



—y/3bi — \/3a —bi — a 6i + a y/Sbi + \/3a 

\/3ai + \/3b —ai — b —ai — b y/3ai + \/3b 

ai — b —\/3{bii + ai) \/6{bii + ai) b — ai 



(G5a) 



(G5b) 



(G5c) 
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